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ON THE ORIGIN OF THE IRREGULAR SATELLITES 
By GERARD P. KUIPER 
YERKES OBSERVATORY, UNIVERSITY OF CHICAGO 
Communicated September 25, 1951 


Thirty satellites are known in the solar system. They fall into three 
classes: 


1. The regular satellites. 
2. The irregular satellites. 
3. The moon. 


The regular satellites are the two of Mars, the inner five of Jupiter, the 
inner seven of Saturn and the five of Uranus, 19 in all. The regular 
satellites have nearly circular orbits, their motion is direct (in the same 
sense as the planetary rotation) and the inclination with respect to the 
planetary equators are all less than 2°. Furthermore, the spacings of these 
satellites are roughly in a geometrical progression, as is true for the planets 
around the sun. More accurately, the spacings appear to depend on the 
masses of the satellites in essentially the same manner as is true for the 


planetary system; i.e., the systems of regular satellites are homologs of the 
planetary system.' This fact has led? to an interpretation of the origin of 
both the planetary system and of the regular satellites in terms of tidally 
stable proto-planets and proto-satellites, formed in each case from a disk- 


shaped nebula by the action of gravitational instability. 

The moon is an exceptional object. Its large mass, 1/81 of its primary, 
indicates that it is not an ordinary satellite. For all other satellites, and for 
the planets to the sun, the mass ratio is less than 10-*. The lunar composi- 
tion (density of olivine, 3.3; absence of an iron core) further indicates that 
the moon was formed as a (win planet with the earth. 

The ten irregular satellites are listed in table 1. The second column gives 
the semimajor axis of the satellite orbit a, in terms of the action radius of 
the planet, R,. The latter is given by the interpolation formula’ 


log (R4/ap) = 0.318 log wp — 0.327, 10-' > wp > 10-5. (1) 


Here a, is the semimajor axis of the planet orbit and w= Mp (M+ Mp), 
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the fractional mass ratio. For Jupiter, R4/ap is 0.0516, for Saturn 0.0352 
and for Neptune 0.0204. R, is the largest distance from the planet to which 
a satellite may go without being lost to the planet, provided its velocity in 
the coordinate frame rotating with the planet is zero at this maximum dis- 
tance. For finite velocities, the maximum distance is less. It is seen that 
the stability of J VIII is precarious, its apojove being at the distance 
0.80 Ra. 

The eccentricity e is exceptionally large for Nereid, while Triton (close to 
the planet) has a circular orbit. The orbital inclinations, 7, are given with 
respect to the planet’s equator except where indicated otherwise. The 
angle of Jupiter's equator being only 3° with respect to its orbit, the six 
Jupiter satellites have moderately large (20-30°) inclinations with respect to 
either reference system. Triton’s present inclination with respect to Nep- 
tune’s orbit is 40°, but this value has little physical meaning since it varies 
(from 40° to 0°) owing to precession of the satellite orbit on the planetary 
equator. The last column shows the direction of motion: D = direct, R = 
retrograde. Itis seen that there are five satellites in each class. 


TABLE 1 
NAMI as/RA ¢ 
JVI 29 0.16 
JVI 29 
JX 29 
J VIL 5M 
J 1X 59 
J XI 56 
lapetus 07 
Phoebe 26 
Triton OO4 
Nereid 0.06 


“With respect to planet's orbital plane 


There might, perhaps, be some doubt as to whether Iapetus should be 
included among the irregular satellites. It was classed among them because 
all clearly regular satellites have 1 < 2°. Quite possibly Hyperion, with 
e = 0.104, is an irregular rather than a regular satellite in spite of its small 7 
value (<1°). All clearly regular satellites have e < 0.02. The strong per- 
turbations which Hyperion experiences from Titan makes the original 


nature of its orbit uncertain. 

Elsewhere‘ the writer has proposed two alternative explanations for the 
retrograde satellites: (1) it was found that collisions between the outer parts 
of consecutive proto-planets can cause retrograde motion of the detached 
parts with respect to one of the two colliding proto-planets; (2) the de- 
crease of mass on the part of all developing proto-planets will cause the loss 
of certain satellites formed before the planetary mass reached its ultimate 
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minimum value. The writer wishes now to withdraw hypothesis (1), as in- 
effective, and put forward the second hypothesis as an explanation of all 
irregular satellites, retrograde and direct. 

The mechanism proposed operates as follows. Let the planet decrease 
its mass by the factor D after a given satellite is formed. It will be seen 
below that D may be as large as 10 or even 100. The law of areas will not be 
affected by this decrease of mass (since it merely depends on the force being 
central); but the energy of the system, and hence ds, will be. The areal 
velocity is ra%/P = constant, if P is the period; while Kepler's third law 
gives ai/P? = cMs/D, if Mg is the mass of the planet at the time of satellite 
formation. Elimination of P gives as/D = constant. The ratio as/R, 
will increase because of two effects: a, ~ D and Ry, by equation (1), is 
proportional to D~°-*, Therefore 


as/R, = constant: D!':*, (2) 


The value of D is limited by D,,,,, being the total mass decrease beginning 
with the early proto-planet stage. D,,,, 1s 10° for the Earth and Venus, 
about 20 for Jupiter, 90 for Saturn and 120 for Neptune.® D values as 
large as 10-10? may therefore oceur for satellites formed rather early in the 
proto-planet history. Such satellites will be lost unless they formed at 
very small ay/R, values, i.e., very close to the center of the proto-planet. 
Viscosity of the proto-planet will change the dynamical argument only for 
satellites embedded in its denser parts; they will move out more slowly 
than computed by equation (2). 

A satellite that has thus been shed by a parent planet will continue to 
move around the sun in an orbit closely resembling that of the planet. It is 
expected to be confined approximately to the zone ap + Ky. It is im- 
probable that the planet just reached its final (present) mass when the 
satellite left it; the general case will be one in which the planet continues to 
lose mass, i.e., one in which its capture cross-section was still large. Sooner 
or later the lost satellite may collide with the proto-planet and be captured 
by it. Such capture may result either in direct or retrograde motion around 


the planet, depending on the geometry of the collision. A collision leading 
to retrograde motion would offer somewhat more resistance to the body 
than one leading to direct motion, so that among the recaptured satellites 


some preference for retrograde orbits is expected. The numbers, 5 and 5, 
are too small to define the empirical ratio with precision. 

Because of the flatness of the proto-planet, no very high inclinations of 
the recaptured satellites are expected. The two groups of Jupiter satellites 
might actually represent only two captures, with each mass having been split 
into three parts, not unlike the break-up of meteorites upon entering the 
earth's atmosphere. 

The eccentricities of the recaptured satellites are expected to have been 
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initially high since as seen from the planet these bodies entered with slightly 
hyperbolic orbits, transformed into elliptical ones only by friction. Re- 
peated passages through the denser parts of the proto-planets would there- 
upon have rounded off the orbits. In one case, Triton, the penetration has 
evidently continued to deeper and deeper layers, causing a very small value 
of ag/R, and a practically circular orbit. Since a large reduction in eccen- 
tricity implies a large increase in mass by accretion, the comparatively large 
mass of Triton is consistent with its low orbital eccentricity. Only the sub- 
sequent dissolution of the proto-planet must have prevented Triton (mov- 
ing against the rotation of the gases) from falling on the planet proper. 
After the satellites were recaptured, their new ds values will have in- 


creased again. In principle, the process of loss and recapture could have 


occurred more than once. When the planetary masses finally stabilized, the 
recaptured satellites might have been at almost any value of as/R, between 


Qand 1. This is actually the case (table 1). 

It is sometimes assumed that the so-called Trojan asteroids were formed 
in the asteroid ring (2 <a < 4 astr. units) and moved to their present a 
values (5.2) by perturbations of Jupiter. This assumption is contrary to 
the conclusion based on Jacobi’s integral for the system Sun-Jupiter- 
Asteroid, namely that all asteroids with a < 4.2 and with small orbital ec- 
centricities (e) and inclinations (7) must remain there.* For finite e and i 
values the stability limit is somewhat less than 4.2 astr. units. 

Only fourteen asteroids are known with a > 4.0: (1) Thule, a = 4.25, 
e = 0.06 andi = 2° (with respect to the ecliptic) is evidently a border-line 
case, though at present somewhat outside the stability limit; (2) twelve 
Trojans, a = 5.1-5.2; and (3) Hidalgo, a = 5.79, e = 0.65, 1 = 43°. It is 
suggested that the Trojans and Ilidalgo are shedded satellites of Jupiter. 
Hidalgo can come close to Jupiter; this fact, and the large e and 7 values, 
support the view that Hidalgo’s orbit resulted from a close approach to 
Jupiter. It is of interest, in this connection, that the absolute magnitudes of 
the six irregular satellites of Jupiter, of the twelve Trojans and of Hidalgo 
are comparable; they correspond to diameters of the order of 50 km. 

If the Trojans are in fact shedded satellites, they are expected to be nearly 
spherical and show little variation in light owing to rotation (< 0.1-0.2 
mag.). If, instead, they are typical asteroids, they are expected to show light 
variations of about 0.4 mag. amplitude and with periods between 5 and 10 
hrs.; such results were found empirically for a sample of the asteroids.’ It 
is hoped that photometric tests now scheduled will decide between these al- 
ternatives. 

' Kuiper, G. P., Astrophys. J., 109, 308 (1949). 

2 Kuiper, G. P., these PROCEEDINGS, 37, 1 and 383 (1951). 

3 This equation is equivalent to, but somewhat more accurate than, equation (17) 
on p. 385 of this Volume. A full discussion of its derivation will be given elsewhere. 
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The formula is precise for w = 10°? and 10° and leads to errors of 1% in Ry at p = 
10-!, 10-3 and 107§, 

* Astrophysics, edited by J. A. Hynek, McGraw-Hill Book Co., Inc., New York, N. Y., 
1951, p. 405. 

§ These PROCEEDINGS, 37, 386, table 1, and 389 (1951) 

® Eg. C.L. Charlier, Die Mechanik des Himmels, Vol. 2, W. de Gruyter, Berlin, 1927, 
Sec. Ed., p. 293 

? Harris, D. L., and Kuiper, G. P., Astrophys. J. (in press). 


THK ORIGIN OF PHOSPHORUS IN 15 BACTERIOPHAGE OF 
ESCHERICHIA COLI 


By Louts W. LABAwW 


LABORATORY OF PHysicAL BroLoGy, NATIONAL INSTITUTE OF ARTHRITIS AND META 
BOLIC DISEASES, NATIONAL INSTITUTES OF HEALTH, U. S. Pustic HEALTH SERVICE, 
FEDERAL SECURITY AGENCY, BETHESDA 14, MARYLAND 


Communicated by R. W. G. Wyckoff, September 28, 1951 


The use of isotopes has given considerable information concerning the 
origin of chemical components of Escherichia colt bacteriophages.'~° Our 
recent studies employing radioactive phosphorus have shown that, though 
various percentages of the phosphorus they contain originate in the host 
bacterium before infection, all seven bacteriophages of the T system obtain 
from the host substantially the same absolute amount of phosphorus per 
infectious unit. In an effort to interpret this equality a more detailed 
study of the origin of the phosphorus in T5 bacteriophage has now been 
carried out. These results show that the bacterial phosphorus incorporated 
into T5 bacteriophage has been taken up by the host, /:. coli, within the two 
minutes before infection. The mean generation time of /. coli in the me- 
dium used is 23 minutes and the burst time for T5 bacteriophage is 40 
minutes. These time relations indicate that, contrary to conclusions drawn 
from work with other T bacteriophages, the bacterial phosphorus going into 
the T5 bacteriophage does not come from the desoxyribonucleic acid of the 
host bacteria. 

Methods.-The F. coli cultures were taken from the same B strain used in 
previous studies with radioactive phosphorus; they were maintained as 
described in previous papers.” * The T5 bacteriophage was that used in a 


previous investigation® and was originally obtained from Professor M. H. 
Adams. It was found necessary to supplement the tryptose broth with 
CaCl to give the 10~*.7 concentration needed for a satisfactory yield of 
T5.” The concentration of radioactivity in broth labeled with P*? was ad- 
justed to give approximately 10° counts per minute per ml. and the data 


were all normalized to this figure. 
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The T5 bacteriophage was purified by sedimenting once from the filtered 
lysates and twice more from 0.85% NaCl. The three ultracentrifugations, 
at an average of 12,000 g. for 30 minutes in the 1.58-cm. diameter tubes of a 
No. 40 Spinco angle head, were interspersed with low-speed clearings in an 
angle centrifuge. 

The number of infectious T5 units was determined by counting plaques 
on horse meat agar plates inoculated with mixtures of /. coli and appro- 
priate dilutions of T5 after incubating the plates for 1S hrs. at 37°C. For 
these titrations the diluted T5 suspensions were allowed to stand with /. 
coli for 15 minutes to facilitate adsorption before the mixture was spread 
on horse meat agar and the plates dried in a 43°C. oven. Surviving bac- 
teria which presumably had failed to adsorb particles of the bacteriophage 
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FIGURE 1 





The amount of labeled phosphorus incorporated into the T5 bacteriophage of 
F. coli as a function of the time of addition of P*? to the culture. These experi- 
mental points are averages where more than one value was measured (see figure 2). 


were determined from dilution platings of the /. coli-T5 mixture on tryp- 
tose phosphate agar. This tryptose phosphate agar does not contain 
enough calcium to support TS multiplication and hence an accurate measure 
of the uninfected bacteria can be obtained if the plates are dried rapidly. 
As usual the radioactive assay was made by measuring dried pellets of 
bacteriophage under an end-window counter. By making these pellets 
from an aliquot fraction of the final T5 suspension, their radioactivity is 
correlated with the number of infectious particles they contain. 
Results.--In order to get information about the state of the phosphorus 
that T5 bacteriophage inherits from its host, the inorganic phosphorus in 
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the culture medium was labeled at different times before and after addition 
of infecting phage to the culture of growing bacteria. Small volumes (130 
ml.) of the parent /. colt, or E. coli + T5, cultures were labeled at various 
times. All the bacterial cultures had bacteriophage added in the ratio of 10 
to 15 per bacterium when the bacterial count had reached 7 X 10’ per ml. 
Lysis was complete in two hours with an average burst size of about 200. 
The T5 bacteriophage in each 130-ml. volume was then purified. The num- 
ber of counts per minute per infectious unit was calculated from the radio- 
active assay and two independent plaque count assays involving the final 
T5 pellet. This was plotted as a function of the time difference between the 
introduction of P*? and the addition of T5. As shown in figure | averages 
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An expansion of part of the curve in figure 1, showing all the experimental points in 
the region —30 to 60 minutes. The two small graphs show how the absorption rate 
affects the results. The data from each separate experiment are plotted with a different 


symbol 


have been plotted where more than one experimental point was obtained 
for a given time difference; the entire plot is a composite of points from five 
experiments of twelve points each. Figure 2 1s a similar curve containing all 
points derived from experiments in which P** was introduced between 30 
minutes before the inoculation with bacteriophage to 60 minutes there 
after. 

These figures show that the amount of radioactive phosphorus incor- 
porated into the T5 bacteriophage is independent of the time of labeling 
the bacteria up to about two minutes before infection with T5., 
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The radiophosphorus in the bacteria, on the other hand, is dependent 
upon the time of labeling since the bacterial phosphorus doubles every 23 
minutes during the part of the logarithmic growth phase used in these ex- 
periments.’ Thus at the time of infection, about 90% of the /. coli phos- 
phorus had been derived from the medium during the previous 85 minutes, 
whereas only 13% had come from the medium in the 5-minute interval be- 
fore infection. The lack of dependence of the bacterial phosphorus in T5 
bacteriophage on the percentage of bacterial phosphorus labeled at the 
time of infection indicates that this phosphorus can scatcely have been 
derived from the desoxyribonucleic acid of the bacteria. 

The right-hand parts of the curves in figures | and 2, where the labeling 
has been introduced after infection, show first a relatively sharp drop to the 
level given in the earlier paper® as the uptake after infection by T5 of in- 
organic phosphorus from the medium. Following this for the first nine 
minutes after infection, there is practically no additional phosphorus uptake 
from the medium; but this level state 1s followed by a further drop continu- 
ing with reduced slope throughout the two hours of the experiment. The 
low uptake of radioactivity after the end of the burst time attests to the lack 
of appreciable adsorption of P** on the bacteriophage, especially since the 
bacteriophage purified for these last experimental points remained in the 
radioactive labeled lysates for over 24 hrs. 

The sharp drop and subsequent leveling off of the curve in the interval 
from 0 to 10 minutes after addition of bacteriophage seems to indicate that 
adsorption of T5 blocks any further utilization for bacteriophage produc- 
tion of phosphorus the bacteria can take up during this period. The £. 
coli-T5 complex begins to take up phosphorus from the medium for bac- 
teriophage production after about 10 minutes, and this continues until 
lysis occurs. The role of T5 adsorption in arresting uptake by the bacteria 
is shown by the difference between the curves in the two small diagrams of 
figure 2. Here the experiment involving slower adsorption gave higher phos- 
phorus values for P** additions in the interval from zero to ten minutes;' 
this would be expected if the bacteria take up phosphorus normally until in- 
fected, 

The ordinates of figures | and 2 may be approximately converted to 
absolute phosphorus contents by using the ratio of (chemical) inorganic 
phosphorus to radioactive phosphorus per ml. of culture medium. In 
these experiments this was measured as: 

104yP/ml. yP 


: ; = 1.04 X 1075 ns 
10° counts, min. / ml. count/muin, 


Hence the phosphorus content per T5 bacteriophage would be about 3.9 X 
10°" g. in agreement with the previous measurement.® 
Discussion. — It is tempting to use the new experimental data to speculate 
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on the nature of the bacterial phosphorus contributed to bacteriophage. 
The present experiment shows that the bacterial phosphorus reaching T5 
is taken from the medium no more than two minutes before infection. It 
therefore can scarcely be the phosphorus in the bacterial nucleic acids. This 
is substantiated by the following calculation. The bacterial phosphorus 
contributed to T5 is equivalent to all of the normal phosphorus uptake by 
uninfected bacteria for a period of about 1'/, minutes. This figure is ob- 
tained by dividing the product of bacterial phosphorus per T5 bacterio 
phage, 1.2 *K 10°" g. P, T5, and the average burst size of 200, by the up- 
take of phosphorus per minute by similarly cultured /. coli {1.6 &K 10>" g. 
P/bacterium/min.].* If then one computes the normal bacterial DNA pro- 
duction per minute taking the DNA phosphorus as 13°) of the total bac- 
terial phosphorus"! one finds that this requires an uptake of only 2 K 10~'¢ 
g. P per minute. This, for two minutes, would only be 17°} of the bacterio- 
phage phosphorus that is contributed by the bacteria. 

Acknowledgment. This work was carried out in the laboratory of Dr. 
Ralph W. G. Wyckoff whom the author wishes to thank for many helpful 
discussions. 

The radioactive phosphorus used in this investigation was supplied by the 
Oak Ridge National Laboratory on allocation from the Isotopes Division, 
United States Atomic Energy Commission. 

Summary. The amount of labeled phosphorus incorporated into the T5 
bacteriophage of /scherichia coli has been measured as a function of the 
time of addition of radioactive phosphorus to the medium both before and 
after infection. The results show that the part of the T5 bacteriophage 
phosphorus which comes from the bacteria before the addition of the virus 
to the broth culture is taken from the medium in the two minutes preceding 
infection. This precludes the bacterial desoxyribonucleic acid as the major 
source of the bacterial phosphorus contributed to T5 bacteriophage. 
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THE EXTENDED ONSAGER EQUATION 

By Epwarp G. BAKER AND CHARLES A. KRAUS 

METCALF RESEARCH LABORATORY, BROWN UNIVERSITY 
Communicated September 18, 1951 

The extended equation of Onsager! may be written: 
Ao’ — Ao = BC + DC log C, 
where 
Ao’ = (A + BYC)/(l — aVC), 2) 


B and D are empirical constants and a and @ are the usual constants of the 
simple Onsager equation. If both B and D are zero, values of Ao’ as calcu- 
lated from the experimental data will be equal to A, within the limit of ex- 
perimental error. The equation thus expresses the deviation from the 
simple Onsager equation as a function of concentration. If D is zero and B 
is positive, we have the equation of Shedlovsky® which holds at lower con- 
centrations for a considerable number of salts. If ion association occurs, the 
deviation curve is still linear but B becomes negative if association is 
marked. 

For a considerable number of salts, particularly of potassium, investiga- 
tors have obtained a better fit to the experimental data by employing equa- 
tion (1) as written.* However, up to now, measurements have not been 
carried out with sufficient precision and analyzed in such manner as to de- 
lineate the form of the deviation curve nor the precision with which the 


equation may be fitted to the experimental data. 

The plot of equation (1) passes through a minimum and is tangent to the 
axis of ordinates at C = 0. To find the condition for the minimum we set 
d( Ao’ — Ao)/dC = 0; solving the resulting equation and converting to 


common logarithms, we obtain: 
0.43 + log C, = —B/D (3) 


where C,, is the concentration at the minimum. Above the minimum, there 
is a concentration, C*, where BC and DC log C are equal and of opposite 
sign. At this concentration, Ay’ = Ao. We thus have 


log C* = —B/D. (4) 
Combining (3) and (4) we have 


0.43 + log Cy, = log C* 


C* = 2.7 C, 
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If the concentration of the minimum can be determined, it is thus possible 
to evaluate Ay without graphical extrapolation since Ay’ = A» at concentra- 
tion C*. 

The value of D is most readily obtained by graphical means. Let us 
transfer the constant Ay to the right-hand side of equation (1) and let us 
write B/D = n. Then equation (1) becomes 


Ao’ = Ao + D(n + log C)C. (6) 
Making use of equation (4), equation (6) may be written in the simple and 
convenient form: 
Ao’ = Ao + DC log c Cc (7) 
If we plot values of Ao’ against values of C log C/C*, we should obtain a 
straight line whose slope is equal to D and whose intercept with the ordi- 
nate corresponding to log C/C* = 0 yields the value of Ap. 
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Plot for the evaluation of D and dy». 


We present an analysis of our results with potassium nitrate at concen 
trations below 5 K 10°*N. From a plot of Ay’ vs. C, Cn was found to be 
14 X 10°-'N. This yields C* = 37.8 K 10°-'N. At this concentration, Ag’ 
(from a Ao’ vs. C plot) is 145.12 and this is, therefore, the value of Ay. For 
B/D =n = —log C*, we obtain the value of 2.4225. 

With this value of n, a plot was constructed of Ao’ against C log C/C* 
(equation (7)). This plot is shown in figure |. The points are numbered in 
order from the most dilute to the most concentrated solution. All points but 
the first one lie on a straight line with a deviation less than 0.005%. This 
plot is very sensitive; the departure of the first point from the straight line 
corresponds to an experimental error of only 0.028%. For Ao we obtain 
145.110 as against 145.120 as found above. The lower value is the more re- 
liable since it was obtained by means of a linear interpolation. 

From the plot we obtain the values: D = 89.7, Ayo = 145.110, B = nD = 
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217.3. With the aid of these constants, we have the numerical conductance 
equation: 


Ao’ = 145.110 + 217.3C + 89.7C log C (8) 


Using this equation, we have computed values of Ao’ at the measured con- 
centrations. The results are shown in the following table. 


POTASSIUM NITRATE AT 25° 


POINT 
NUMBER Cc x 10 Ao’ (exp.) Ae’ (cale.) AAo’ 


$9432 145.116 145. O81 +-0 035 
844448 145.071 145 067 +0). 004 
76103 145.056 145.060 —0.004 
6750 145.052 145.057 —0.005 
1318 145.071 145.067 +0. 004 
(272 145.158 145.155 +0.003 





A,-A,: 217. 3C + 89.7C logC 
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Plot of equation (1). 


In figure 2, the smooth curve represents the numerical equation; the 
circles represent experimental values of Ao’; their radii correspond to an 
error of 0.01 A-unit in the experimental conductance values. With the ex- 
ception of the first point, the deviation of the experimental values from the 
equation is less than 0.0050). It would appear, therefore, that for potas- 
sium nitrate, the extended Onsager equation reproduces the observed con- 
ductanices very precisely. 

The deviation curves of all potassium salts measured (chloride, bromide, 
iodide, nitrate and bromate) exhibit minima, as also do cesium picrate and 
chloride and sodium picrate. The plot for sodium iodide is linear with posi- 
tive slope. Quaternary ammonium salts of large anions yield linear plots 
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with positive slope but their salts with small anions yield linear plots with 
negative slope. 

' Onsager, L., Physik. Z., 28, 277 (1927); Onsager, L., and Fuoss, R. M., J. Phys. 
Chem., 36, 2689 (1932); Fuoss, R. M., Physik. Z., 35, 59 (1935). 

? Shedlovsky, T., J. Am. Chem. Soc., 54, 405 (1932). 

* Shedlovsky, T., Brown, A. S:, and MaclInnes, D. A., Trans. Electrochem. Soc., 66, 
165 (1934); Benson, G. C., and Gordon, A. R., J. Chem. Phys., 13, 473 (1948); Monk, 
C. B., J. Am. Chem. Soc., 70, 328 (1948); Owens, B. B., Jbid., 61, 1393 (1939) 


CONFIGURATIONS OF POLYPEPTIDE CHAINS WITH FAVORED 
ORIENTATIONS AROUND SINGLE BONDS: TWO NEW PLEATED 
SHEETS 


By Linus PAULING AND ROBERT B. CoREY 


GATES AND CRELLIN LABORATORIES OF CHEMISTRY,* CALIFORNIA INSTITUTE OF 
TECHNOLOGY, PASADENA, CALIFORNIA 


Communicated September 4, 1951 


In recent papers we have described several configurations of polypeptide 
chains with interatomic distances, bond angles, and other structural fea- 
tures as indicated by the studies in these Laboratories of the structure of 
crystals of amino acids, simple peptides, and related substances, and have 
presented evidence for their presence in synthetic polypeptides, fibrous 
proteins, and globular proteins.'~* The requirements that we have imposed 
for a satisfactory polypeptide configuration, in addition to the correct bond 
distances and bond angles, are that each amide group be planar, with 
either the cis configuration or the trans configuration about the C’—N bond 
(which has nearly 50 per cent double-bond character), and that each car- 
bonyl and imino group (except for proline or hydroxyproline residues) be 
involved in the formation of a hydrogen bond with N--H---O distance 
approximately 2.8 A and with the oxygen atom nearly on the N--H axis. 
The only other structura] parameters involved in a configuration of poly- 
peptide chains are the orientations around the NC and C— C’ single bonds. 
In the following paragraphs we discuss the question of the relative sta- 
bility of structures with different values of these orientational parameters. 

We are interested in the potential function for orientation around a single 
bond between the @ carbon atom, which forms three single bonds in addi- 
tion to the bond under consideration, and either the nitrogen atom, which 
forms a single bond to its hydrogen atom and a bond with somewhat less 
than 50 per cent double-bond character to the amide carbon atom, or the 
carbon atom C’, which forms a bond with oxygen which has somewhat more 
than 50 per cent double-bond character and a bond with nitrogen with 
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somewhat less than 50 per cent double-bond character. We shall discuss 
only polypeptides containing residues with the L configuration. The po- 
tential function would be expected to have six minima, approximately 60° 
apart, and six maxima. If the three single bonds on the a carbon atom were 
equivalent and the two bonds formed by the other atoms were equivalent, 
as for example in nitromethane, H;C--NO», the potential function would 
have the form shown in figure | (a), with equivalent minima and equivalent 
maxima. For propylene, Hy;C--CH-—-CH,, the potential function has either 
the form shown in figure | (4) or that in figure 1 (c); either the minima are 
divided into non-equivalent sets of three each, or the maxima are so 
divided. Curve Jb corresponds to the situation in which the stable con- 
figurations are those with one of the bonds of the @ carbon atom coplanar 


yet Vt wr ea NT 


III 
JAI V/ \ I\/ /\ 


Bred WI WIEN 


360 
FIGURE | 


Forms of the potential function for 
orientation around a single bond con 
necting a tetrahedral carbon atom to a 
planar group 


with the CH--CHz, group (the two sets representing the configuration cis 
to the double bond and trans to the double bond, respectively) and curve 
Ic represents the situation in which the stable configurations are those with 
one of the single bonds of the a carbon atom in a plane normal to the plane 
of the CH=-CHp group. 

The height of the potential hump for nitromethane has been independ- 
ently determined by three investigators, by measurements of heat capacity 
and of entropy,'®~' their average value being 0.6 keal. mole ~' (Pitzer and 
Gwinn, 0.8 keal. mole~'; DeVries and Collins, 0.6 keal. mole~'; Jones and 
Giauque, 0.5 kcal. mole~'). For propylene the height of the potential 
hump is larger; the value 0.8 kcal. mole~' or less originally suggested by 
Pitzer" has been replaced by the value 2.1 keal. mole~', which is supported 
by several other investigators.''~'’ In neither case have the thermodynamic 
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data provided a choice between the alternatives represented by curves 1b 


and Ic. For the bonds that we are considering all of the maxima and all of 
the minima would be different, and the potential function would be expected 
to have a form such as that shown by curve 1d. 

It seems reasonable to expect that the average energy differences of the 
maxima and minima for an amide would be somewhere between 0.6 kcal. 
mole~' and 2.1 keal. mole~', perhaps about 1.0 kcal. mole~', with consider- 
ably smaller differences among the minima themselves. This energy quan- 
tity, | keal. mole~', is sufficiently large to cause configurations represented 
by the minima in the curve to be significantly favored over those repre- 
sented by the maxima. On the other hand, in crystals of amino acids and 
related substances in which hydrogen bonds, with energy as great as 8 kcal. 
mole~', are formed the metrical requirements might be such that the most 
stable structure, involving the best conditions for hydrogen-bond forma- 


FIGURE 2 


A diagrammatic representation of two 
possible stable orientations of the bonds 
of the a-carbon atom with respect to the 
plane of the amide group. The aC atom 


" 


is directly beneath C’; an adjacent single 


bond (tapered lines) of the aC atom lies 
either (a) in the plane of the amide group 
or (6) in the plane perpendicular to it 


tion, is one in which the molecule is distorted by rotation around the single 
bonds to the a carbon atom into a configuration of instability with respect 
to the orientational potential functions. With proteins in an aqueous en- 
vironment the effective energy of hydrogen bonds is not so great, inasmuch 
as the difference between the energy of the system with N—H.- --O hydro- 
gen bonds surrounded by water and a system with the N--H group and the 
O atom forming hydrogen bonds with water molecules may be no more than 
about 2 kcal. mole.~' Hence the favored orientations around the bonds to 
the a carbon atoms may be of more frequent occurrence in proteins than in 
simple amides. 

The first question to be answered is whether the adjacent single bond on 
the a carbon atom lies in the plane of the amide group or in the plane per- 
pendicular to this plane (Fig. 2). The evidence from the known structures 
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of simple substances is not conclusive, but indicates that the stable con- 
figuration is the coplanar one. In N-acetylglycine'* the bond from the a 
carbon atom to the carboxyl carbon is very nearly coplanar with the amide 
group, the carboxyl carbon lying cis to the N--H group. In #-glycylgly- 
cine!’ the same configuration is shown by the carboxyl carbon atom and the 
amide group. Inaddition,ineach of these twosubstances the carboxyl group 
itself lies in the same plane as the amide group (to within 0.03 A in 8-glycyl- 
glycine), so that the single bond from the amide nitrogen atom to the a 
carbon atom has assumed a coplanar orientation with respect to the car- 
boxyl group. In the benzyl penicillinate ion” the bond between the carbon 
atom of the benzene ring and the @ carbon atom adjacent to the carbonyl 
carbon atom of the amide group is coplanar with the amide group, and 
trans to the CO group. There is also a carbon-carbon single bond in- 
volving the @ carbon atom attached to the nitrogen atom of the amide 
group that is coplanar with the amide group, and cis to the N--H group; 
the significance of this orientation is decreased, however, by the fact that 
this carbon-carbon bond is part of a four-membered ring. In 8-glycylgly- 
cine, on the other hand, the amino nitrogen atom is rotated by 27° out of 
the plane of the amide group. We consider it likely that this corresponds to 
a strained configuration, and that the larger amount of evidence indicating 
the stability of the orientation corresponding to a coplanar arrangement of 
one of the single bonds of the a carbon atom and the amide group is to be 
accepted. 

Further evidence is provided by the structures of synthetic polypeptides 
and proteins. The most common configuration of polypeptide chains 
seems to be that of the 3.7-residue helix, which we have proposed as repre- 
senting the structure of synthetic polypeptides,‘ hair, horn, quill, muscle, 
and related fibrous proteins,’ and hemoglobin;® valuable confirmation of 
the assignment of this structure to the substances has been provided by 
Perutz’s discovery that the 1.5 A reflection corresponding to the fiber-axis 
length per residue, which had been observed for porcupine quill, can also be 
obtained from the plane normal to the helical axis for the other substances.*! 
In this helix the C--H bond and the C—R bond (that is, the bond to the ~ 
carbon atom of the R group) lie within a few degrees of the plane of an ad- 
jacent amide group. In the three-chain structure that we have proposed for 
collagen and gelatin’ five of the six C--H and C-—R bonds lie very near 
the planar configuration, the other one being rotated through nearly 30°, 
presumably into the unfavored configuration to which it is constrained by 
the nature of the structure. 

Accepting the conclusion indicated by the preponderance of the evidence 
that the minima for the rotational potential functions for the two bonds 
connecting the a carbon atom with adjacent atoms in the polypeptide chain 


correspond to coplanar configurations, we give an exhaustive discussion in 





Vou. 37, 1951 CHEMISTRY: PAULING AND CORE} 


the following paragraphs of configurations of this type in which all of the 
amide groups are of the trans type and are equivalent. 

There are six favored orientations about the C—-N bond, which we may 
describe as CH trans to NH, CR trans to NH, CC’ trans to NH, CH cis to 
NH, CR cis to NH, and CC’ cis to NH. Similarly, we may describe the six 
favored orientations around the C—C’ bond as CH trans to CO, CR trans to 
CO, CN trans to CO, CH cis to CO, CR cis to CO, and CN cis to CO. There 
are accordingly 36 different configurations of the polypeptide chains of L- 
amino acid residues corresponding to these favored orientations; all 36 
configurations are different. They are tabulated below in four groups: 


a, 


[act c ‘C ypc ac q 
= / | 


FIGURE 3 


A diagrammatic representation of a hydrogen-bonded 
layer structure of polypeptide chains with alternate 
chains oppositely oriented, showing steric hindrance be- 
tween 8-carbon atoms of adjacent chains 


trans-trans, cis-cis, trans-cis, and cis-trans. It might be anticipated that 
the configurations in which CH, CR, or CC’ (or CN) is trans to the bond 
with the greater amount of double-bond character, of the two bonds once 
removed, would be either more stable or Jess stable than the configurations 
in which one of these bonds is cis to the bond with the greater amount of 
double-bond character. Accordingly we would expect that either the nine 
cis-trans structures or the nine trans-cis structures would be somewhat more 
stable than the other 27 structures; however, the differences in stability 
are presumably small enough to be unimportant in comparison with the 
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TABLE 1] 


CONFIGURATIONS OF EQUIVALENT TRANS AMIDE GROUPS WITH FAVORED ORIENTATIONS 
ABouT C--N anp C--C’ Bonps 


A. Trans-trans Configurations: 

1. CH trans to NH, CH trans to CO. | Pitch 6.68 A, 2 residues per turn. Lateral 
CO and NH suited to forming the two pleated sheets described below. 

2. CR trans to NH, CR trans to CO. Like 1. Suited to forming the two 
pleated sheets, but prevented (except for polyglycine) by steric hindrance 
of R groups 

CC’ trans to NH, CN trans to CO. Ruled out by steric hindrance of NH 
and O 

CH trans to NH, CR trans toCO. Pitch 6.2 A, 4 residues per turn. Lateral- 
diagonal CO and NH, suited to forming a tetragonal framework with 
a = 4.8 A, but prevented by steric hindrance of R groups. 

5. CR trans to NH, CH trans to CO. Like 4. 

CH trans to NH, CN trans toCO. Ruled out by steric hindrance of NH and 
N of adjacent group. 
CR trans to NH, CN trans to CO. Like 6 
CC’ trans to NH, CH trans to CO. Ruled out by steric hindrance of O and O. 
CC’ trans to NH, CR trans to CO. Like 8. 
Cis-cis Configurations: 

10. CH cis to NH, CH cis to CO. Ruled out by steric hindrance of NH and O. 

11. CR cis to NH, CR cis to CO. Like 10 

12. CC’ cis to NH, CN cis to CO. Pitch 7.30 A, 2 residues per turn. Lateral 
CO and NH, suited to forming the planar sheet with alternate chains re 


versed in direction, but prevented (except for polyglycine) by steric hind- 


rance of R groups 

3. CH cis to NH, CR cis to CO. Pitch 5.5 A, 3.7 residues per turn. Axial CO 
and NH, suited to forming the a helix with intramolecular hydrogen bonds. 

14. CR cis to NH, CH cis to CO. Like 13; the a helix. 

15. CH cis to NH, CN cis to CO. Ruled out by steric hindrance of O and O. 

16. CReis to NH, CN cis to CO. Like 15 

17. CC’ cis to NH, CH cis to CO. Pitch 10.4 A, 4 residues per turn. Lateral 
CO and NH, with H of NH near axis, and CO projecting out. Not suited 
to formation of intramolecular or intermolecular hydrogen bonds (except 
possibly for polyglycine only, in a tetragonal framework with a = 4.8 A) 

18. CC’ cis to NH, CR cis to CO. Like 17. 

Trans-cis Configurations: 

19. CH transto NH, CH cisto CO. Pitch 5.2 A, 4.8 residues per turn. Diagonal 
CO and NH, with NH directed toward axis, CO outward Not suited to 
formation of hydrogen-bonded sheet or framework 

20. CR trans to NH, CR cis to CO. Like 19 

21. CC’ trans to NH, CN cis to CO. Ruled out by steric hindrance of O and O. 

22. CH trans to NH, CR cis to CO. Pitch 8.0 A, 2.8 residues per turn. Not 
suited to formation of hydrogen-bonded sheet or framework. 

283. CR trans to NH, CH cis to CO. Like 22. 

24. CH trans to NH, CN cis to CO. Pitch 9.6 A, 2.7 residues per turn. Lateral 
CO and NH, not suited to formation of hydrogen-bonded sheet or frame- 
work 

25. CR trans to NH, CN cis to CO. Like 24. 
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26. CC’ trans to NH, CH cis to CO. Ruled out by steric hindrance of O and 
C’ of adjacent amide group 
27. CC’ trans to NH, CReisto CO. Like 26. 
Cis-trans Configurations: 
28. CH cis to NH, CH trans to CO. Pitch 5.2 A, 4.8 residues per turn (similar 
to 19). Not suited to formation of hydrogen-bonded sheet or framework 
29. CR cis to NH, CR trans to CO Like 28. 
30. CC’ cis to NH, CN trans to CO. Ruled out by steric hindrance of NH and 
NH 
CH cis to NH, CR trans to CO. Pitch 7.6 A, 2.7 residues per turn. Lateral 
CO and NH, not suited to formation of hydrogen-bonded sheet or frame- 
work. 
CR cis to NH, CH trans to CO Like 31 
CH cis to NH, CN trans to CO. Ruled out by steric hindrance of NH and 
adjacent amide group 
CR cis to NH, CN trans to CO. Like 33 
CC’ cis to NH, CH trans to CO. Pitch 9.5 A, 2.7 residues per turn. Lateral 
CO and NH. Not suited to formation of hydrogen-bonded sheet or frame- 
work 
36. CC’ cis to NH, CR trans to CO. Like 35. 


differences due to other structural features. (Note that the orientation 
trans to CO is trans to the bond with the larger amount of double-bond 
character, whereas the orientation cis to NH 1s trans to the bond with the 


larger amount of double-bond character.) 
In table | we give for each of the 36 structures, all of which are helical, 


the pitch of the helix (axial length per turn), the number of residues per 
turn, the approximate orientation of the CO and NH groups (relative to 
axial, diagonal, or lateral), and the suitability of the struc 


the helical axis 
ture to forming hydrogen bonds. Hydrogen bonds may be intramolecular 
or intermolecular. In the latter case either a hydrogen-bonded sheet or a 
hydrogen-bonded three-dimensional framework may be formed. A sheet 
can be formed only when the helix has a twofold screw axis (or a onefold 
screw axis, which does not occur), and a three-dimensional framework, 
under our requirements of equivalence of all of the amide groups, only 
when the helix has a threefold, fourfold, or sixfold serew axis. In de- 
scribing the configurations we have allowed deviations of a few degrees 
from the favored orientations around the single bonds, in order to permit 
the formation of the suitable hydrogen-bonded structures. 

From table | we see that the only configurations of interest are 1, 12, 13, 
and 14. The last two configurations represent the @ helix, with two dif 
ferent sets of positions of the side chains. Configuration 12, a completely 
extended polypeptide chain, with identity distance 7.30 A, is suited to the 
formation of a planar sheet, as shown in figure 4, except that steric hindrance 
between adjacent chains prevents the formation of this sheet in case that 
the side chain R is anything but a hydrogen — that is, this planar sheet could 


be formed only by polyglycine. 
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The remaining configuration, 1, is of great interest —it seems likely that it 
is present in many proteins with the #-keratin configuration, The pre- 
dicted identity distance along the chain axis is 6.68 A, corresponding to 3.34 
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FIGURE 4 


A diagrammatic representation of the antiparallel-chain 
pleated sheet structure. 


\ 


\ 
HCBC HCBC HCBC 


FIGURE 5 
A diagrammatic representation of the  parallel-chain 
pleated sheet structure 
A per residue. This agrees very well with the residue length reported by 
Astbury and Street for proteins with the 6-keratin structure, 3.32 A. 
We have previously suggested for 6 keratin a pleated-sheet configuration 
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of polypeptide chains in which each residue is converted into the following 
one by the operation of a glide plane of symmetry. With this configuration 
the orientation around the bonds N-—C and C—C’ is not that which we con- 
sider to be stable—the stable orientations lead to a residue length in the 
chain of only 2.84 A, which is not found experimentally for any protein of 
the 8-keratin type. We now think that it is unlikely that this pleated sheet 
is an important configuration for proteins. 

We have found that polypeptide chains with the configuration | of Table 
1 and residue length 3.34 A can be arranged to form sheets with lateral 
hydrogen bonds in either one of two ways. In the first way alternate chains 
in the sheet have opposite directions. The sheet is similar to the first 
pleated sheet that we described in having @ carbon atoms at the same level 


TABLE 2 


Atomic CooRDINATES FOR THE ANTIPARALLEL-CHAIN PLEATED SHERT AND PARALLEL- 
CHAIN PLEATED SHEET 


Antiparallel-Chain Pleated Sheet: 
Atomic coordinates: x, y, 2; %,'/2 + y, 2; '/2 , §, 8; + x, "/2 — y, 2. 
ao = 946A 
by = 6.68 A (direction of chain axis) 
c = 1A (arbitrary) 
Parallel-Chain Pleated Sheet: 
Atomic Coordinates: x’, y, 2; & 
a = 4.73 A 
bo = O68 A 
o =1A 


x y x’ 
0 009 0.000 OF O18 
—0.086 0.184 25 O72 
0.0538 0 316 25 106 
0.183 Q). 293 16 366 
0.044 0.000 -2.39 O88 


in adjacent chains displaced to the front and the rear of the plane of the 
sheet, so that it is appropriate to describe it as a pleated sheet. We sug- 
gest that this sheet be called the antiparallel-chain pleated sheet, the 
original one being called perhaps the polar pleated sheet, inasmuch as all of 
the carbonyl groups are oriented in the same direction, either to the right or 
to the left, the imino groups being oppositely directed. 

The second way of arranging the chains of type 1 into a sheet involves 
parallel chains, which also are arranged to form a pleated sheet. Diagram- 
matic representations of the antiparallel-chain pleated sheet and parallel- 


chain pleated sheet are shown in figures 4 and 5, and drawings of these 


structures in figures 6 and 7, respectively 
Atomic coordinates for the two pleated sheets are given in table 2. The 
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antiparallel-chain pleated sheet, considered as made of L-amino acid resi- 
dues, has orthorhombic symmetry, its symmetry elements being a twofold 
screw axis in the direction of the chains, a twofold screw axis normal to this 
direction and in the plane of the sheet, and a twofold axis normal to the 
plane of the sheet. The parallel-chain pleated sheet has only a twofold 
screw axis in the direction of the chains. 

In the parallel-chain pleated sheet the oxygen atom lies nearly in the 
N--H direction, The deviation from the N—-H axis would be zero in case 
that the bond angle C--N--H were equal to 125°. It is likely that this 
bond angle is about 116°, and that there is accordingly a bend in the 
hydrogen bond of about 10°. This is not so great as to lead to significant in- 
stability 
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FIGURE 6 


Drawing representing the anti-parallel-chain pleated sheet structure. 


It is interesting to note that the ability of polypeptide chains with con- 
figuration | to form lateral hydrogen bonds with either a parallel or an anti- 
paralle] adjacent chain strongly suggests that proteins with the 8-keratin 
structure may exist in which there is a randomness in sequence of parallel 
and antiparallel chains. We plan to discuss this question later, in connec- 
tion with a comparison of observed and calculated intensities of x-ray re- 
flections for synthetic polypeptides and proteins of the 6-keratin type. 

It is to be noted that our 5.1-residue helix' “* is not represented in table 1, 
although it satisfies the conditions of containing trans amide groups, and 
of having all amide groups equivalent. The reason for its failure to appear 
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as a result of the considerations of the present paper is that it represents 
unfavorable orientation about both the N-—C bond and the C--C’ bond; 
the orientation about each bond is such as to correspond to a rotation of 
about 30° for both C-—-H and C—R from the position of coplanarity with 
each of the two amide groups adjacent to the a carbon atom. If we accept 
the assumption expressed above that these orientations represent an in- 
stability of about | kcal. mole ~' relative to the favored orientations, then 
the 5.1-residue helix would be predicted to be less stable than the 3.7-residue 
helix by 2 kcal. mole~! per residue. There is no indication from the struc- 
tures that a corresponding extra stability would be conferred on the 5.1- 
residue helix by, for example, the presence of N—-H- --O hydrogen bonds of 
abnormally great stability, the absence of steric hindrance between side 
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FIGURE 7 


Drawing representing the parallel-chain pleated sheet structure, 


chains of adjacent turns of the helix which might be present in the 3.7- 
residue configuration and make it less stable, or any other structural fea- 
ture —indeed, as was pointed out in an earlier paper,” the cylindrical cavity 
down the center of the 5.1-residue helix would be expected to lead to an 
extra instability for this configuration by decreasing the intramolecular van 
der Waals attraction. We thus are led to the conclusion that the 5.1- 
residue helix is a less stable configuration of polypeptide chains than the 
3.7-residue helix, and that it would not be expected to occur very often as 
an important feature of the structure of proteins. In particular, we are 
provided with grounds for withdrawing our tentative suggestion’ that pro- 
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teins of the a-keratin type undergo transformation to the 5.1 residue con- 
figuration on supercontraction. This suggestion was made because the 
amount of contraction involved in the change from the 3.7-residue helix 


to the 5.1-residue helix, about 35 per cent, is just that reported to occur in 
the process of supercontraction. We now think that it is likely that, as 
suggested by Astbury, the process of supercentraction involves merely a 


disorientation of a-keratin molecules, without significant change in the 
configuration of the individual polypeptide chains. 

This investigation was aided by grants from The Rockefeller Foundation, 
The National Foundation for Infantile Paralysis, and The National Insti- 
tutes of Health, Public Health Service. 
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HYBRIDS BETWEEN LYCOPERSICON ESCULENTUM MILL. 
AND SOLANUM LYCOPERSICOIDES DUN. 


By CHarLes M. Rick* 
DIVISION OF TRUCK Crops, UNIVERSITY OF CALIFORNIA 
Communicated by Karl Sax, September 14, 1951 


The garden tomate (Lycopersicon esculentum Mill.) has been successfully 
hybridized with all other species of Lycopersicon that have been grown for 
experimental purposes.'? Its crossability relations with the large parent 
genus, Solanum, on the other hand, are largely unknown, although a 
strong barrier to crossing with certain species might be inferred from the 
following observations. Solanum melongena, S. muricatum, S. nigrum, S. 
tuberosum, and other species, especially of the tropics, are often grown close 
to tomatoes, or are found as weeds in tomato fields. Despite this oppor 
tunity for crossing, spontaneous hybrids between the tomato and_ these 
species have apparently never appeared, and apparently no such hybrid 
has been produced artificially. The present example of a tomato-night- 
shade hybrid was encountered in an exploratory survey of the relations be- 
tween the two genera. 

The tomato parent of the hybrids is Pearson, the leading variety of the 
California tomato industry. It is typical of the standard or normal-type L. 
esculentum in the genetic sense except that it is homozygous for the sp gene, 
which conditions determinate plant habit. 

The nightshade parent is Solanum lycopersicoides Dun., a perennial 
native to the south coast of Peru. It has been classed by Bitter*® in the 
section 7Tuberarium, the group of Solanum species usually considered to be 
most closely allied to the tomatoes.4- Even among species of this section, S. 
lycopersicoides shows a unique resemblance, as its name implies, to the 
tomatoes in such characteristics as its yellow flowers, nature of leaf dissec- 
tion, and lack of tubers. 

Taxonomically S. /ycopersicoides is a true Solanum in having the following 
features that distinguish it clearly from Lycopersicon: its anthers lack sterile 
tips, dehisce terminally and do not join laterally to each other to form an 
anther tube. 

Plants of S. lycopersicoides were acquired? as collection No. 30382 of the 
University of California Botanical Garden Third Expedition to the Andes, 
1942. Seeds had been collected near Candarave, Peru (Dept. Tacna, Prov. 
Tarata). Plants in cultures at Davis resemble very closely the herbarium 


specimen of the original collection (Spec. No. 695,139 of the University of 


California Herbarium). 
The F; hybrids between these two species were obtained with considerably 
greater ease than those between L. esculentum and L. peruvianum (L.) Mill. 
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The parents were crossed in both directions, but since fruits are set er- 
ratically on the nightshade parent even after selfing and intraspecific 
crossing, the tomato was mostly used as the pistillate parent. Following 
pollinations of the latter kind, small fruits are set containing ovules that 
develop to an intermediate size. A few small seeds reach maturity. 
Since these seeds did not germinate, 30-day fruits were opened under sterile 
conditions, the ovules dissected, and the endosperm-embryo contents 
placed in sterile cultures of 0.89 agar supplemented with Hoagland and 
Snyder’s solution,® minor elements, and 2% dextrose. In these cultures 
embryos grow slowly at first, many being deformed, but after an initial 
period of slow growth, they develop very rapidly. 

Fifteen embryos were successfully cultured in this fashion, and eight of 
these have been transferred to soil. Further vegetative growth in the green- 
house and field is not attended with any delay or abnormality; in fact, 
under field conditions the hybrids grow with much greater vigor than either 
parent. 

The two parent species differ in a great number of readily distinguished 
characters. Comparisons between the parents and F, hybrids for 25 of these 
traits have been made. Each sample consisted of 20 or more measurements. 
Only three characters of L. esculentum-. day-length response, fragrance of 
flower and relative calyx segment length— are dominant in the hybrids. Of 
the remaining 22, five of S. lycopersicoides are completely dominant in the 
hybrids, the other 17 being intermediate but usually showing greater re- 
semblance: to S. lycopersicoides. In respect to the key characters of the 
anther the hybrid would very likely be annotated as a Solanum; in respect 
to general appearance the hybrid seems like a magnified S. lycopersicoides. 
Considering the morphology, the subsequently described sterility, and 
chromosome behavior of these plants, their hybrid origin can scarcely be 
doubted. 

The fertility of the hybrids is extremely low. Of the pollen grains formed, 
only 0.5°% are stained by acetocarmine, and many of these appear other- 
wise defective. Ovule fertility as assayed in self-pollinations and back- 
crosses to both parents, though not yet completely analyzed, is at least as 
low as the pollen fertility. If L. esculentum is used as the pistillate parent in 
backcrosses, occasional ovules develop, and preliminary embryo cultures 
have yielded one seedling. It is suspected that with perseverance a back- 
cross progeny might be grown, and that characters might be transferred 
from SS. lycopersicoides to L. esculentum. 

Chromosome behavior can explain at least part of the very Jow fertility of 
the F, hybrids. Beth parents possess 12 pairs of chromosomes, those of SS. 
lycopersicoides appearing larger at metaphsase than those of L. esculentum. 
The cytology of two hybrid seedlings was studied, each revealing the ex- 
pected 24 chromosomes. In 79 figures at metaphase of the first meiotic 
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division a mean of 14.8 of the 24 chromosomes was paired, mostly as bi- 
valents, only a few trivalents and quadrivalents being observed. The num- 
ber of paired chromosomes varied from 8 to 24. The unpaired chromo- 
somes lag and usually fail to enter telophase nuclei. A mean of 9.6 chias- 
mata was observed in 30 pollen-mother cells at diakinesis. Occasional 
figures exhibit chromosome bridges. 

Considering this cytological picture and the nearly perfect pairing of 
chromosomes in all hybrids between species of Lycopersicon, the chromo- 
somes of L. esculentum seem to be closer in homology to those of other 
tomato species than to those of S. /ycopersicoides. 

The relatively low degree of chromosome pairing in the hybrid suggests 
that some degree of fertility might be restored in amphidiploids. This ex- 
pectation has been confirmed to the extent that tetraploids obtained from 
colchicine treatment exhibit 67%) pollen fertility and nearly complete pair- 
ing of chromosomes, mostly as bivalents. 

The feasibility of obtaining Lycopersicon-Solanum hybrids opens a new 
and interesting opportunity for further research in the genetics and breed- 
ing of the tomato. In addition to the many known features of the tomato 
that favor it for cytogenetic research, a genetic link with the Solanum 


species is now available that provides an avenue for exchange of genes be- 
tween the two genera. Too little is known of the traits of S. /ycopersicoides 
to permit any evaluation of it as a parent for practical breeding purposes; 


in fact, it is not even known whether its fruits contain any of the poisonous 
alkaloids that are often encountered in fruits of other nightshade species. 
At least, it is known to be immune to, or tolerant of, tobacco mosaic, for, 
despite frequent contact with infected tomato plants during two years of 
clonal culture, no symptoms have appeared. 

It might be appropriate to consider here the bearing of certain facts about 
the hybrids upon the relations between the parent genera, although an ex- 
tensive consideration of phylogeny here would be premature. A reasonably 
close relationship between the parents of the present hybrids must be as 
sumed because the latter were obtained with relative ease and because they 
have great vegetative vigor. Although the hybrids are highly sterile, 
backcrosses to the tomato parent and possibly to both parents can suc- 
ceed, thereby providing a medium of exchange of genes between the two 
parents. Such features are encountered in many known hybrids between 
species of the same genus -for example, Nicotiana. In fact, most inter- 
specific crosses within plant genera fail, and many that succeed yield hybrids 
of low viability and of even greater sterility than the present examples. 
The latter thus behave genetically more like hybrids between species of the 
same genus rather than of different genera, and one is led to question the 
separation of Lycopersicon and Solanum into different genera. 

Summary,—F, hybrids have been obtained between Lycopersicon 
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esculentum and Solanum lycopersicoides by resort to embryo culture. The 
hybrids are more vigorous vegetatively than either parent and in the 
majority of morphological characters resemble the nightshade parent 
more closely than the tomato parent. Pollen and ovule fertility are very 
low, but it is possible to secure progeny from the backcross to the tomato 


parent. The hybrids possess 24 chromosomes, the same number present in 
both parents. Slightly more than half of the chromosomes pair in the dip- 
loid hybrids, but in colchicine-induced amphidiploids nearly all of the 
chromosomes pair, chiefly as bivalents and pollen fertility is greatly in- 
creased. 

* This study was supported in part by a grant from the John Simon Guggenheim 
Memorial Foundation. 
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A MAJORANT METHOD FOR NON-LINEAR PARTIAL 
DIFFERENTIAL EQUATIONS 


By S. BERGMAN AND M. SCHIFFER* 


MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
AND STANFORD UNIVERSITY 


Communicated by T. Y. Thomas, September 18, 1951 


1. Let ) bea finite simply connected domain in the (x, y)-plane bounded 
by a smooth curve C. Consider the non-linear partial differential equation 
fake 


Au = P(x, yu + ¥ a,(x, y)u’, = ? (1) 
? Ov" 


where P(x, y) is a positive continuous function in D, and the coefficients 
a,(x, y) are arbitrary continuous functions in D satisfying the inequalities 

l bad cai , 

|a,(x,y)| $< K-» K > 0,/ > 0 constants. 

v! 
This last inequality guarantees the convergence of the series on the right- 
hand side of (1) for arbitrary values of u. 
We prescribe the boundary values y(s) of u(x, y) on C(O €sS Ls= 

length parameter on C) and ask for a solution of (1) assuming these bound- 
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ary valueson C. The value u(x, y) of the solution at a given point (x, y) ¢ D 
is a functional of the boundary values u(s) and it is natural to try to de- 
velop the solution u(x, y) into a functional power series of u(s).2 For sake 
of simplicity, we reduce this development immediately to an ordinary power 
series, as follows. Let u(x, y; ¢) be that solution of (1) which has on C the 
boundary values /u(s) and let us set up: 


u(x, y; tl) = Zz u,(x, ye (3) 
v=l 
We assume the following boundary values for the coefficients u,(x, y) of 
this development: 


u(x, vy) = m(s), u(x, y) = O for vy > 1, (x, y) e C. (4) 


2. We show at first that if a development (3) exists at all it is uniquely 
determined by the differential equation and the boundary conditions (4). 
In fact, introducing (3) into (1) and comparing equal powers of ¢, we obtain 


Au, = P(x, ym, u, = p(s) on C, 
Au, = P(x, y)u2 + a(x, yur’, te = OonC, 


(5) 


Au, = P(x, y)tt. + Eq(Ge, .. . Gay Miy+ +> Met), ty = Oon C. 


Here /:, is a well-defined polynomial, linear, homogeneous in the a,’s and of 
degree n in the u,’s with positive coefficients which is easily computed for 
any given positive integer n. 

The determination of the coefficients u,(x, y) is thus reduced to the well- 
known theory of the linear partial differential equation :'~* 


Au = P(x, y)u, P(x, y) > Oin D. (6) 


Let G(x, y; £, 7) be the Green's function for this equation with respect to the 
domain D, We may then solve the inhomogeneous linear differential equa- 
tion for u,(m > 1) in the form: 


u(x, 9) = — S SoGlx, y; &, n)En(ae(, 0), . .  Unalé, n))dé dn. (7) 


Since “,(x, y) is uniquely determined and can be easily expressed in terms of 
Green's function also: 


. ae 
u(x, vy) = Je uls) 5 G(x, y; &, n) ds, vy = interior normal, (8) 
Vv 


we can successively compute all coefficients of the development (3). The 
treatment of the linear partial differential equation (6) may be reduced to 
the theory of orthogonal particular solutions and to the theory of the kernel 
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function. From a numerical point of view, the development (3) together 
with formulas (7) and (8) seems, therefore, to be quite convenient. 

3. The decisive problem in our investigation is now the convergence of 
the series (3). We will study it by means of a majorant procedure based 
upon the dependence of the Green's function G(x, y; & 7) upon the co- 
efficient P(x, y) and upon the domain D. We remark at first that Green's 


: ae ish ¥e F OG ‘ 
function is positive in D and that, therefore, > >Oon C. Let 
Vv 


max |mu(s)| = A. (9) 
Oss SL j 
We replace our previous boundary value problem by the following: To find 


a solution for the partial differential equation 


A (Tho 
Au* = P(x, yju* — K >> ' (10) 


2 wv! 


which has on C the constant boundary value A. Let u*(x, y; t) = 


> u, (x, yt’ be the solution of (10) with the boundary value tA on C, and 
yeu 
let kK; have an analogous definition to /,. In view of (8), the coefficient 
u(x, y) in the new problem will be positive and at each point larger than 
the absolute value of the previous ™(x, y). Since each new — E;, will be 
larger than the absolute value of the corresponding /, and since Green's 
function is positive, we will have 
un(x, y) 2 | un(x, y) |. (11) 
The series (3) will therefore, converge absolutely wherever the corresponding 
series for u*(x, y; t) does converge. We may, therefore, restrict ourselves to 
the convergence proof for this new series. 
Let 
min P(x, y) = m. (12) 
(x,y) «D 
It is well known that Green's function for a fixed domain D decreases every- 
where if the coefficient P(x, y) of the equation (6) is replaced by one which 
is larger everywhere in D.’ Hence, let G,,(x, y; &, 7) be the Green’s function 
of D with respect to the partial differential equation 


Au = mu; (13) 
we have then the inequality 
G(x, y; é, n) ~ Ga lx, J; £, ), 


and since both Green's functions vanish on C, we also have 
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G . Ga 


nic. 5 
Ov Ov ” (45) 


Let us replace the partial differential equation (10) by the simpler equa 
tion 
_ @, (lv)’” 
Av = mv —- K >> a (16) 


9 


v < 


with the same boundary condition v = A on C. Defining in analogy to (3) 
v(x, y;t) = D> o,(x, ye’, (17) 
v=l 


we easily conclude from (7), (8), (14) and (15) that 


‘ 
v(x, y) u,(x, y) 2 | U(X, y) |. (1S) 


v 


Let, finally, 'y be a circle of radius R which contains the domain D. 
Consider the radially symmetric solution of the linear partial differential 
equation 


Aw = mw, w = w(r), (19) 


which is normalized in such a way that its minimum on C is exactly one. 
Let \(C) be the value of w(R); this number depends on the geometry of C 
and can easily be expressed in terms of Bessel functions. 

We now pose the following boundary value problem: To find a solution 
L(x, y) of the partial differential equation (16) in I’, which has on the cir- 
cumference the constant boundary value AA(C). We define again the 
power series 


U(x, yt) = ¥ U(x, yt’. (20) 
v= 


We remark that U(x, y) is defined by the equation: 

AU, = ml, U,; = AX(C) on I'g. 
Comparing this equation with (19), we find 

U(x, y) = Aw(r) 

and since we have 

Av; = my v(x, y) = A < Ux, y) on C, 
we derive by the maximum principle 

v(x, vy) S U(x, y) in D. (24) 


We denote Green's function for the circular domain and the differential 
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equation (13) by g(x, y; & 9). Since Green’s function increases mono- 


tonically with the domain, we have the inequality 
g(x, ¥; & 1) Z Gm(x, ¥, & 0), (x, vy) € D. (25) 


Since the 1’, as well as the v, are obtained by recursion formulas (7) based 
on the Green’s functions, we have by (24) and (25): 
’ | ! . 
U(x, y) 2 vfx,.y) 2 | u(x, y) |. (26) 
Thus, the convergence of the series (3) will be ensured, if we can prove 
the convergence of the corresponding series in the particular boundary 
value problem for U(x, y). But in this case, U(x, y) will only depend on the 
radius r to the center of the circle !'y. We may reduce the whole question 
to a problem in ordinary differential equations. In fact, U(x, y) = U(r) 
satisfies the ordinary differential equation 
poet =e cg ww 
(r) + r) =mU — Klev — 1 —1U] 27 
dr? r dr ail 
with the boundary conditions U(O) = finite, U(R) = AX(C). It is easily 
seen that the solution of (27) depends analytically upon U(R) if this value 
is small enough. Hence, the series development (3) will converge for 


| > rr. . 
max |u| = A sufliciently small. The existence theorem for the boundary 
value problem with respect to the non-linear partial differential equation 
(1) has been reduced to the study of a boundary value problem for the or- 


dinary differential equation (27). 

It should be pointed out that in the theory of partial differential equa- 
tions in a given domain D the majorant method is based on a comparison of 
the Green's functions for different domains and different differential equa- 
tions. It is a procedure in the large in contradistinction to the Cauchy- 
Kowalevsky method of majorant power series which guarantees solutions 
only locally and is, therefore, not applicable in boundary value problems. 

The method of majorants can readily be extended to other problems in 
non-linear partial differential equations of elliptic type. It may be applied 
to the case of more than two independent variables and may also be ex- 
tended to the case of boundary value problems of the second and the third 
type. In the latter problems, the Neumann's and the Robin's functions will 
play analogous roles to the Green’s function in the above development. 

In order to study the breakdown of the series development (3) and the 
uniqueness of the solutions for given boundary values, one has to study the 
analytic function (3) in dependence of its argument ¢ in the complex plane 
and to continue it analytically beyond its singular points on the real axis, 


if possi ble. 
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ON FINITE GROUPS WITH TWO INDEPENDENT GENERATORS. 
WI. EXPONENTIAL SUBSTITUTIONS 
By JESSE DOUGLAS 
CocumsBia University, New York Crry 


Communicated September 18, 1951 


The two previous notes of this series have appeared in earlier issues of 


these PROCEEDINGS.' This one is a direct continuation. 
16. Conjugates of a Given Linear Substitution. —Let 


g(x) = px (mod n), (B,n)=1 — (16.1) 


be a given linear substitution,’ not reducing to the identity: 8 # 1 (mod 
n).® We deal now with the problem of finding all its conjugates 6. 

The range m of @ must be a multiple of the period s of ¢, as defined by: 
Bp = 1 (mod n), where s has its minimum positive value. We write, as 
usual, m = su. wis one of the unknowns in the determination of 6; it will 
be the modulus in most of the congruences occurring hereafter. 

Another unknown is the 6-parameter g: 


g O(s)/s (mod un) (16.2) 


We know that g must be prime to u; let A denote the exponent to which g 
belongs modulo u, i.e., the minimum positive exponent for which the con- 
gruence 

g” (mod 1) 

is valid. g has the properties 


t 


g” l, ee es g” l (mod wm) 
(see Theorem XV, Note IIT). The last relation is obviously implied for all ¢ 


by its validity for / E: 
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| (mod u). (16.3) 
If « = (6 — 1, nm), we may summarize by saying that 
g=!l (mod u), (16.4) 


or that A must divide e. 
The conditions for conjugate special substitutions 6, ¢ are: 


6(0) = 0 (mod m), ¢(0) = 0 (mod n), 
6"(x) =x (mod m), ¢y”" (x) =x (mod n), (16.5) 


7,0 = 0°”, Ti¢” = ¢*”. (16.6) 


(Cf. Note I, $3.) 

Now g”"(x) = B"x = B“x = x (mod n);_ thus the second condition 
(16.5) is met. 

T\¢” is the transformation 


x g(x+1)—- ¢(0) =P (x+ 1) - P= Ppx (mod n), 


aly e 
"” the transformation 


and ¢ 
x’ = Bx (mod n). 
The second condition (16.6) is therefore equivalent to 
ly) = y (mod s). (16.6a) 
Accordingly, we may write 


O(x) =x + sP(x) (mod m), (16.7) 


¥(x) being an unknown function, whose determination is precisely our goal. 
Of course, we must have 


y(0) =0 (mod u) (16.8) 


in order that 6(0) = O (mod m). 
A further necessary condition is that @ admit the translation 7); this is 
equivalent to 


V(x + s) = V(x) + Ws) (mod 14). 
From (16.2), (16.7), 
=1+ ys) (mod u); 
hence 
W(x + 5s) = W(x) + g — 1 (mod wu). (16.10) 


By induction on ./, we infer 
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v(x + Ms) = P(x) + Mg — 1) (mod xu). (16.11) 


in order that 6 really represent a substitution on the range m, we must first 
require that if x = y (mod m), then @(x) = @(y) (mod m); equivalent to this 
(by (16.7)) is ¥(x) = Wy) (mod u), which follows from (16.11) when we re- 
place M by Mu. Conversely, 6(x) = 6(y) (mod m) must imply x = y 
(mod m). Now, in fact, 6(x) = 6(y) (mod m) implies @(x) = @(y) (mod 
s), which (by (16.6a)) implies x = y (mod s),orx = y + Ms. Therefore 
(by (16.7)): y + Ms + sly + Ms) = y + s¥(y) (mod m), M + p(y + 
Ms) = ¥(y) (mod u), Mg = 0 (mod u) (by (16.11)), M 0 (mod u) 
since g is prime to u, finally x = y (mod m), Q.E.D. 

It remains only to render valid the first conditions in (16.5) and (16.6) by 
proper choice of the function y. 

We have 


67(x) =x + s¥(x) + s¥(x + sP(x)) 
=x + s¥(x) + stWl(x) + (g — I(x)} 
=x + s(x) (1 + g) (mod m), 


and by induction on y we readily establish: 


. 
(x) == x + s(x) _ teak (mod m). (16.12) 
oo] 

If g = 1, the fraction (g’ — 1)/(g — 1), being =1+¢+e¢?+...+ 2's 
to be interpreted as equal to y. This will be a standing convention through- 
out, whatever letters appear in place of g and y; e.g., in formula (16.18), the 
fraction = xifb = 1. 

Using (16.12), we find the first condition (16.6) to be expressed by the 
formula 


l . 
| ix +1) — yl)! =vx (mod m) 


(16.13) 


for all x, y. Now 


By { 


(g” — 1)/(g — 1) = fe” — D/(e — DEA’ — D/A -— VD}, 
and 


I+ e+ pP+... 
El - ge Pr... % (mod u) 
= (g” — 1)/(g -— 1), (16.14) 


where we use the fact that g’ = g (mod u), by (16.3). 


Hence (16.13) becomes 
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I 
| 


y> - | 
tye +1) — wy} =* 


g L gd 


¥(x) (mod w), 


which will hold for all y if and only if it holds for y = 1: 


8 
y — | 
¥(x + 1) — (1) : 4) (mod 1), (16.15) 


4 
Defining 


y(1), 


(mod 1), (16.16) 


we write (16.15) in the form 


y(x + 1) a + by(x) (mod a), (16.17) 
This gives 


v(x + 1) — (x) bi v(x) — Wx — 1)} (mod 1), 


Ay(x) bAy(x — 1), where A(x) = P(x + 1) — W(x). 
Consequently, AY(x) = ab* (x 2 0). Replacing x by y and summing from 
y Otoy x 1, with account taken of (16.8), we get 


x 


v(x) F (mod u) (x = 0). (16.18) 
ae 


(For b = 1, the fraction is to be interpreted as equal to x.) 
Formula (16.9) imposes the condition 


be — 1 


a sia” g- (mod 1), (16.19) 
and with use of the fact that (6 — 1)(¢ — 1) = 0 (mod u) (by (16.16) and 
(16.3)), we easily see that this secures the requirement (16.10). Formula 
(16.11) may be used to extend the definition of ¥(x) to negative values of x. 

It only remains now to take account of the first condition (16.5); by 
(16.12), (16.18), this leads to 


(mod u), 


which is valid for all x if valid for x 1, so giving us our final condition 


(mod ). (16.20) 
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Since g* = 1 (mod u) and X divides both 8 — | and n (ef. 16.4)), we can 
transform (16.16) and (16.20) as follows: 
B-1g° - 

rn g 


(mod 1) 


n - J 


1 (mod 1). (16.22) 
AXg-1 


( 


This is done by a procedure like that which gave us (16.14). 

In summary, to find in the most general way a conjugate substitution 
6 to the linear one ¢ defined in (16.1), we must first select arbitrarily a value 
of u to give us the range m = su of 6, then an integer g, prime to u, which 
obeys (16.4), then define ) by (16.21), then, if possible, determine a so as to 
obey simultaneously the congruences (16.19), (16.22). @ will then be de 
fined by (16.7) with (x) given by (16.18): 

b* 


I 
A(x) v + sa } (mod m). (16.23) 
Te 


By (16.12), we see that the powers of @ are given by the formula 


A(x) x = ; (mod m) (16.28a) 


if g # 1 (mod w), and by the formula 
b* | 
#(x) x + say (mod m) 


’ 
) 


if g 1 (mod 1), 
From this we infer that, if g # 1 (mod «), the period + of 6 is the least 


positive integral exponent such that 
g =! (mod u(g — 1) /(a, u)), (16.238¢ 
while, if g 1 (mod 4), 
u/(a, u). (16.23d) 


r = mis the necessary and sufficient condition that the given ¢ be the 
derivative of 0: g = 4. 
Since by Theorem XVI in Note IT, $9, 
gy | (mod 2), (16.28e) 
one may observe, by noting the definition of following (16.2), that 
Ar. 
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The Case g = 1 (Mod u). — It is convenient to abbreviate the system 
(16.19, 22) to: 
aB=g-— 1, aC = 0 (mod u) (16.24) 


by writing 


» 
_ Be oj 


= (mod u). (16.25) 
a= 


We define 
A = (B, C, 2), D = (C, %). (16.26) 


Obviously, d = (B, D). 
The determination of a is always possible if we select g = 1 (mod w), 
which is in accord with (16.4). 
For the system (16.24) then becomes 


aB =0, aC = (0 (mod u), (16.27) 


where, since \ = | and b = £6 (mod u) (see (16.16) and the convention 
stated just after (16.12)), we have 
B- 1 


B= ; "== n (mod u). (16.28) 
p-1 


The solution for a of (16.27) is: 
a = 0 (mod u/A). (16.29) 


This gives us precisely A values of a that are incongruent modulo uw: 
a= M(u/d), M = 0,1, 2,..., 4 — 1. The substitution 6 (16.23) corre- 
sponding to a general value of \/ is the \/th power of 6’, where 6’ corresponds 
to M = 1. 

We may observe that the period r of is equal to A/(A/, A). Obviously, 
(16.26, 28), Alm, therefore rn, as required by the conditions for conjugacy 
of 0, g. The derivative of 61s: 0,)(x) == Bx (mod r). @is linear in case 6 = | 
(mod r), and is otherwise non-linear, therefore of type 3 (cf. the remarks 
associated with formula (16.34)). 

The Case g # 1 (Mod u). For a general value of g, a brief study of the 
system (16.24) gives the following results. 

Let the integer P be determined by the congruence 

P(B/A) = 1 (mod D/A). (16.30) 
This is always possible, since 4 = (B, D), and therefore B/A, DA are rela- 
tively prime. 

Then 


== P(g — 1)/A (mod u/A), (16.31) 
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provided that (i) P(g — 1)/Ais an integer which, (ii) when substituted for a, 
obeys the congruences 
aB=g— 1, aC = 0 (mod 4). (16.32) 
Otherwise the system of congruences (16,19), (16.22) for @ is inconsistent, 
and there is no substitution @ corresponding to the values selected for u and 
g. In the case of consistency, we have A values of a in arithmetic progres- 
sion: 
a = P(g —1)/A+ M(u/d), M=0,1,2,...4-~1. (16.33) 
For the first and second differences of the function @(x) in (16.23), we 
have 
Ab(x) = 1+ sab’, A’O(x) = sab* (b — 1) (mod m = su). 


It follows that the criterion for the reduction of @(x) to a linear substitution: 
6(x) = (1 + sa)x (mod m), is 


a(b — 1) =0 (mod 4). (16.34) 


Provided that a(b — 1) # 0 (mod wu), the substitution defined in (16.23) 
will be called exponential. Since then @ 1s a conjugate of linear, while not 
itself linear, it is clear that exponential substitutions are precisely those of 


type 3. (Cf. Theorem X, Note I and Theorem XXIII, Note IT.) 
We may also remark that the condition that two exponential substitu- 
tions — (16.23) and the same with a’, b’ in place of a, b--be identical is: 


, 


a=a', ab = a'b’ (mod u), (16.35) 


so that exponential substitutions which look different may really be the 
same. E.g., for all x, 


x + 4(2* — 1) Sx 4+ 2 (1-(—1)*) (mod }2), 
x+1— (-ly=x*4+3 - 1] (mod 8). 


The (¢ + 1)th difference of @(x), g any positive integer, is 
At *'6(x) = sab™(b 1)" (mod m = su); 
hence if there is a positive integer g such that 
a(b 1)’ =0 (mod w), (16.36) 


the function (x) is congruent mod m to a polynomial of degree S q, and 
exactly of degree q if g 1s the minimum positive integral exponent for which 
(16.36) is valid. E.g., the type 3 substitution @ = (0) (2) (4) (6) (1357) is 
represented exponentially by the formula 


O(x) =x +1 —-— (-1)* (mod 8), (16.37) 
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subsumed under the general formula (16.23) by writing s = 2,u = 4,a = 1, 
b 1. Since here 


a(b — 1)? = 4 0 (mod 4), 


whereas a (b — 1) 2 #0 (mod u), (x) must really be quadratic; in 
fact one easily verifies: 
OX) mat Xf 2x* (mod 8). (16.38) 
17. I:xamples.—(a) we consider the problem of finding all the conjugates 
of 
¢g(x) -x (mod n), (372) 
the range n being arbitrary. 
We have 8B = —1,s = 2,€ = (2, m). 
If mis odd, we must have e = 1; therefore g= 1 (modu). This gives us 
b B -1,B=0, C=n(modu). A= (B,C,u) = (n,u). a = 0 (mod 
u/A). In summary, if m is odd, all the conjugates of (17.1) are given by 
the formula 


A(x) x + M (u/A) (1 (—1)*) (mod m = 2u)_ (17.2) 


where \J = 0,1, 2,..., 4 — 1. These are the powers of the particular ex- 


ample 6’ corresponding to JJ = 1. 6’ may be described as the substitution 


which keeps the even numerals fixed and advances each odd one cyclically 
by m/A, so that each cycle of 6’ contains either Aequally spaced odd numerals 
or one even numeral. @ is linear if and only if 247 (u/A) = O (mod w), 
2M = 0 (mod A), = 0 (mod A) (since Ais odd); thus, among the powers 
of 0’ only the identity is linear. If «is prime to n, then A = 1, and (17.2) 
gives only the identity on the range 2u. 

If n is even, « = 2; therefore either g= 1 or g? = 1, g # 1 (modu). By 
choosing g 1, we obtain again the type (17.2). The discussion of the 
case g? == 1, g # 1 can be carried through in each particular instance accord- 
ing to the plan described in the preceding section. E.g., for g(x) = —x 
(mod 4), we find as particular conjugates the substitutions 6, 0~' of example 
(d). 

(b) Conjugates of 


¢o(x) = 3x (mod 8). (17.3) 


on 


We have B = 3,5 = 2,m = 2u,¢€ = (8 — 1, m) = (2,8) = 2, 27?=1 

(mod «). ‘To obtain a particular example, let us take u = 8, g = 3. Then 

m = 16, = 2,6=5, B=6, C= 0 (all mod 8). A = (6, 0,8) = 2, 

D = (0,8) = 8. 3P = 1 (mod 4), P = 3 (mod 4), a = 3-2/2 = 3 (mod 4) 
3,7 (mod 8). Choosing a = 3 (mod 8), we have 
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A(x) x 4 "(3 (mod 16). (17.4) 
Expressed in cycles: 6 = (0) (S) (1, 7, 9, 15) (3, 13, 11, 5) (2, 6) (4, 12) 
(10, 14). The choice a == 7 (mod 8) gives 6-'. We observe that the basic 
valuea = P(g — 1)/d = 3is an integer which obeys Ba g-1,a@@=0 
(mod u), as required. 

The linear conjugates 6 of ¢: A(x) ax (mod m) are precisely those of 
period 2 on an even range: a? 1 (mod m),m = 2u. (Cf. Theorem XXVI, 
Note IT.) 

18. Antiderivatives of a Given Linear Substitution. A conjugate @ of a 


given substitution ¢ does not necessarily have ¢ for its derivative. We can 
only say: ¢(x) = @(x) (mod r), where r, the period of @ and the range of 
6,, is a divisor of n, the range of ¢. As has been stated, in order that 


¢ = 6, r must be equal to n. 

The problem of antiderivatives—given ¢, find 6 such that 6 yg 1s im- 
portant because from a given substitution it enables us to find one of the 
next higher type. 

We consider here the problem of the antiderivatives of a given linear sub- 
stitution (16.1). The following is the procedure leading to the desired re 
sult. 

Let ¢ have the period s, i.e., let 6 belong to the exponent s, modulo n. 
Then m = su. Choose u arbitrarily, then g prime to uw so that g* | 
(mod u), where « = (8 — I,m). Suppose that g belongs to the exponent A, 
modulo u; then Ae, so that An. We calculate: 


yg - l | 
G : G, B : (mod #). 


£ 
(18.1) 


(G, u), 
Then to allow us to proceed further, 
p = n/X must divide (B, d). (18.3) 


Suppose that this condition is met, and let the integer R be defined by the 
congruence 
B‘p R (aod 6). (IN.4) 
Since B is determined mod u, Bp is determined mod ip (d/p)é; hence 


R is indeed determined mod 6, since d/p is an integer by the requirement 


(18.3). 
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Now find /:, the g.c.d. of R and 6; then a requirement for further proce- 


dure is that 
I divide (g — 1)/d (18.5) 


(the latter always an integer). If this condition is fulfilled, let the relatively 
prime integers //, K denote the cofactors of FE in R, 6: 


R = EH, 6= KK. (18.6) 
The congruence for A: 
AHl = (g — 1)/kd (mod A) (18.7) 


will then always have a solution A», unique mod KA, so that the general value 
of A is 


A=A,+ MA (18.8) 


with A/ an arbitrary integer multiplier. 

Now Af must be selected so that A is prime to p:(A, p) = 1. Only the 
values of M in the interval 0 <= MW < pK will come into consideration, and, 
indeed, the values 0 S M < p are decisive for the primality of A to p, since 
Ay + MK and Ag + (AM + p)K are both prime or both not prime relative to 
p. 

In case all this is possible, i.e., if (g — 1)//d is an integer, and if among 
the solutions A of the congruence (18.7) there are values prime to p, then 
the substitution of these values in the formula 


d 
A (mod u) (18.9) 


p 

gives the appropriate values for a. Values of A congruent mod pA’ give 
values of a that are congruent mod dEAK = dé = u; this is the reason that, 
as stated, 7 may be restricted to the values 0 S M < pk. 

The final step is to substitute the values determined for a, 6 in the general 
formula (16.23) for @. 

The case g = | (mod w) is of especial simplicity. The only requirement 
to be met is: 


(Bp — 1)/(p8 1) =O (mod mn). (18.10) 


We must then choose « = a multiple of n, n|u, and have 


u 
B, A- (mod nu) (18.11) 
Nn 


where A is any integer prime to”. Appropriate substitution is then to be 
made in the formula (16.23) for 6. 
example: Antiderivatives of 
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g(x) = Bx (mod 7) (1S.12) 


where 8 — 1, as well as §, is prime ton. This hypothesis obviously insures 
the fulfillment of condition (18.10). 
We have e = I, so that g = 1 (mod uw). Then 


mp — | 
Ox)=x+A (mod m), (18.13) 
np-l 
where the range m of @ must be chosen equal to a multiple of ns, msm. A is 
any integer prime to , which is equal to the period r of @(x). Lf @’(x) 
corresponds to A = I, then (18.13) represents the substitution 6’ to the A 
power, and defines, for (A, m) = 1, all the antiderivatives of ¢. 

When 1 is a prime p, the condition: 8 —— | prime to p, is equivalent to 8 
# | (mod p). With» = p, A an arbitrary non-multiple of p, formula 
(18.13) thus gives us the most general exponential substitution of period p. 

The linear substitutions of period p are defined by 


A(x) = ax (mod m), = l, a # | (mod m). 


For the existence of an a obeying these conditions, it is necessary and suf- 
ficient that p*|m or that m have a prime factor of the form Mp + 1. This 
is seen from the formula for Euler’s g-function: 

m= p*q...r’, g(m) = p’ gr (p l)(g-1)... - 1), 
and from the facet that those residue classes modulo m which consist of inte- 
gers prime to m form, by multiplication, an abelian group of order g(m). 
For an abelian group of order NV contains an element of period equal to a 
specified prime p if and only if p divides N. 

We thus have a census of a// special substitutions of prime period p, since 
such a one must be either linear or exponential (Theorem XXV, Note II). 

The preceding discussion leaves unsettled, for the time being, the ques- 
tion of the existence of an antiderivative for an arbitrarily assigned linear 
substitution——one where the sufficient condition (18.10) is not verified, such 
as: g(x) = 5x (mod 8), or g(x) 4x (mod 9). 

19. Powers of Special Substitutions. We deal with the question: when 
is a given power @* of a special substitution @ of period r also special? The 
answer is: when and only when the rows of index Ay in the array of the 
corresponding primary group [ = (6, 6) form a subgroup of T. This oc 
curs when and only when d 6,(dx) for all x, where d = (A,r). (Cf. Theorem 
XVII, Note II.) A tantamount condition is: d6\(d) for all ¢, i.e., the 
cycle of @, that contains d shall consist entirely of multiples of d. The de 
rivative y = 0,'(x) of 6’ = 6 is then defined by the congruence 


r O,(Ax) ( i 4 
d => d Mo d . 
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always solvable uniquely for the integer y since (1) @,(Ax)/d is an integer, 
(11) A/d is prime to r/d. 

Algorithmically, to apply (19.1), we may divide by d in those cycles of 4, 
which consist of multiples of d (simply omitting the other cycles, in which 
none of the numerals is divisible by d), and then multiply each numeral by 
\’, where (\/d)A’ = 1 (mod r/d), reducing the result to its least non-nega- 
tive residue mod r/d. 

If #6 is an exponential substitution, its derivative 6; is linear: 6,(%) == Bx 
(mod r), and the condition: d,4,(dx) for all x, is obviously fulfilled. We 
have for the derivative of 6’ = 6°:0;(x) Bx (mod r/d). Hence: 

Tureorem XXVIII. Any power 6° of an exponential substitution 0 is a 
special substitution, exponential if 8B # 1 (mod r/d), linear if 8 == 1 (mod 
r/d), where d = (x, r). 

Of course, any power of a linear substitution is linear. But the property 
of the last theorem is not characteristic of exponential and linear substitu- 
tions. E.g., in the case of the following substitution on the range 32: 


6 = (0) (16) (1, 7, 9, 31, 17, 23, 25, 15) (3, 29, 27, 5, 19, 13, 11, 21) 
(2, 6, 18, 22) (10, 30, 26, 14) (4, 12) (8, 24) (20, 28), 


every power of @ is special, but @, having the derivative 
0, = (0) (2) (4) (6) 1 3 5 7%), 


of type 38, is itself of type 4, neither exponential nor linear. 


Ihe necessary and sufficient condition that every power # of a given 


special substitution @ be itself special is that, for every divisor d of r, the 
cycle of 6, containing d shall consist entirely of multiples of d. | This require- 
ment is met, in particular, ifr = p’, p a prime, and p is a fixed numeral of 
6. In the example just given, r = 2%, and the numeral 2 is fixed in 6. 

! Vol. 37, nos. 9, 10 (September, October, 1951) 

2 The parenthesis notation is employed for greatest common divisor. 

§ The contrary case has already been disposed of in Note II, $13. 


MAXIMUM PROPERTIES AND INEQUALITIES FOR THE 
KIGENVALUES OF COMPLETELY CONTINUOUS OPERATORS* 
By Ky Fan 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF NOTRE DAME 
Communicated by John von Neumann, September 8, 1951 

1. Although several extremum properties of the eigenvalues of com- 
pletely continuous (abbreviated to ¢. ¢.) operators in Hilbert space are 
known, nearly all of them deal with a single operator. In this respect, the 
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following theorem of J. von Neumann! is of particular significance, as it in- 
volves two matrices: Let A, B be two square matrices of order n (with complex 
coefficients). If \dij, {xi} are the eigenvalues* of A*A, B*B respectively, then 
n / 

> V dK is the maximum of, tr (UAVB)|, when LU, V run over all unitary 
r= 1 

matrices of ordern. In the present note, Theorem | contains two maximum 
properties (4), (5) for any finite number m of c. c. operators in Hilbert 
space . (4) is a generalization of von Neumann's theorem, while (5) is a 


multiplicative companion of (4). The case m = 1 of (5) includes the im- 
portant inequality of H. Weyl?’ gif... &:}? << Aide... An (@ = 1,2, ...) 
between the eigenvalues | &;{| of ac. c. operator A in Sand the eigenvalues |X| 


of A*A. The case m = 2 of (5) includes the following theorem of A. Horn! 
concerning two c. c. operators A, Bin S: If {ds}, fet, fart are the eigen 
values of A*A, B*B and (AB)*(AB), respectively, then ype... in < Arde 
Ankiky ... Kk, forn = 1,2,.... The proof of our Theorem | is based on 
this result of Horn and the elementary Lemma |. Thus, in the case of 
finite matrices, our proof of (4) differs from von Neumann’s more geometric 
proof, 

For any two c. c. operators A, B in $, inequality (10) in Theorem 2 gives 
a comparison, essentially different from Horn’s, of the eigenvalues of 
(AB)*(AB) with those of A*A and B*B, while inequality (11) compares 
the eigenvalues of (A + B)*(A + B) with those of A*A and B*B. From 
(11) we derive Theorem 3 on asymptotic distribution of the eigenvalues of 
(A + B)*(A + B). In the case of linear integral equations with symmetric 
kernels, both (11) and Theorem 3% are well-known results of Weyl.’ 

As an application of the maximum property (4), we have in Theorem 5 
another inequality comparing the eigenvalues of (A + B)*(A + B) with 
those of A*A, B*B and involving an arbitrary symmetric gauge function. 
This generalizes a result of von Neumann! concerning finite matrices. To 
derive Theorem 5 from (4), we first prove Theorem 4 concerning a general 
inequality for symmetric gauge functions. The geometric meaning of 
Theorem 4 is clear. Analytically it is analogous to a classical theorem of 
R. F. Muirhead,® although the class of functions involved in his theorem is 
entirely different from the class of symmetric gauge functions. 

> Lowmr ta. Beha oa>. eae Oh Ss he... 2 2 
0. If py (lL < 4,7 < mn) are n* non-negative numbers such that 


n 


D py <i forallj, YS py <1 for alli, (1) 
s=1 j= 1 


then 


n 


n n 
‘a } Py ai by < > a, by. 
t=ly=l 1 
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Proof: Letce = YS pybs(1 <i <n). Fora fixed positive integer 
, 


i 
k n 

k<n,letg = py 0<j<n). Thend< q,;< land YS ack. Since 
i, 


b, = bj4, 2 O, we have 


n k-1 k—l k~1 
> why < ¥ qb; + (« -> w) bh. < ¥ (bj — by) + Rd, 
j= | j=l 


j=l j=1 j= 


k k 
Das Dd b, (1<k<n). 
= 1 t=] 


n 


n 
By Abel's lemma, a; > a,4, 2 0 and (3) imply ¥ ac <S > S a,b,. 
1 1 


i=] 
We shall use the above lemma in the following form: 
Lemma In. Let fa,}, |b,{ be two non-increasing infinite sequences of non- 


negative numbers such that > ab <+ ,. If py (i,j = 1, 2,...) arenon- 
i=1 


negative numbers such that >> py < 1 for all j and XY py < 1 for all i, then 
i=) j=l 


> DY pyaihy < ¥ aby. 
sz ly=l i=] 

THEOREM |. Let Ay, As, ..., A,, be m (21) ©. c. operators in Hilbert 
space S. For each] = 1, 2,..., m, let Xj, Ajay. . Ay». « De the eigenvalues” 


. * ‘y *,¢ . 
of A;A,;. Then for any positive integer n: 


Max : (UsAi. .- Unda Be %) | = > tie (4) 
1 1 


n 
Max | Det (U;A,... UpAgxs, Xe) |? = TD Ojde:. ss Xn? (5) 
I<t,ksn t=] 
where, for both maxima, , ls, . . ., Uy, independently run over all unitary 
operators and |x, X», .. ., %,{ runs over all sets of n orthonormal elements in SD. 
Proof: For each j = 1, 2,..., m, let /7, be the non-negative square root of 
A;A;. Then I, is also c. c. and its eigenvalues are | V yy}. Let yy; (1<7 
< n) be n orthonormal elements in such that //,y;; = V diy (lo tcn). 
Let A, = W,//, be the polar factorization of von Neumann.’ Here W, is a 
partially isometric operator and its initial set is the closed linear manifold 
determined by the range of //;. We have then W;A, = H). Since W; is 
also partially isometric and its initial set is the closed linear manifold deter- 
mined by the range of A), we can find a unitary operator |’; in such that 
VjAyyn = WjAyyn = Wyn (1S i< n). Let Z, be a unitary operator 
in ) such that Z)y,, = yj-1,; (I< ¢< n), where we put y,; = vn; Let now 
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U; = Ai V; (1 < 9 < m). Then UA, see U, Aw¥re = (AA se Ant) "Yi 
(1<1< nm). Hence, for the particular choice of U'; = 2,V,(1< 7< m) and 
Xi = Vmi (1S 1< n), the relations 


n 


> (UA, “$y Ee Ui Anti X;) < > (AyAzi sa p ae ‘ (6) 
1 1=1 


Det (UA, wane Ui Aske X,) 2 < Il (AyjAe: ie Ami) (7) 
i=l 


lst, kan 
are satisfied and actually become equalities. 
It remains to show that inequalities (6), (7) hold for any m unitary opera- 


tors U,..., Um and any n orthonormal elements x, ..., x, in . We first 
prove the case m = | of (6): 


n n 

> UiAwn x) |< YS Vw (8) 
= 1 s=1 

Consider the polar factorization of A, = Wh. Let {yu} be a complete 

orthonormal set such that Jy, = WaAn vu (k = 1, 2, 3,...). Then, 

since (U,A\x,, x1) = Uhx, Wil ix,), we have 


p (U\A Xi, X;) = > b V dy (x1, Vie) Vir Wi; Xj) 
i=1 t=lk=l 


and therefore 


n n 


2,0 UA, x) |S DL VM | (Xm yu) |? + 
Lk 1 


by Zz. V dx (Wivu, Uix,) 
i= k i 


t= ' 


By Lemma 1b, each of the last two double sums is< 55 V Ay. Thus (8) is 
i=l 


proved. 
Let now B = U\A,... UnAnm and let {x,{ be the eigenvalues of B*B. 
Since (U,A,)*(U)A,) = AjA), repeated application of Horn’s theorem gives 
n 


< Il (Ay,Az, ee ) or a (n = l 2 eee (9) 
1 


» =» 
i 


nN nN 
which, according to a lemma of Polya,” implies >. VK< Zz. (AgAzi « Ami) 
1 l 1 1 


Hence (6) will result when we apply (8) to B. 
According to Lemmas 2, 3 in an earlier note* by the author, we have for 
any n orthonormal elements x), X2, .. ., Xn/ 


Det (Bx,, x.) |\?<. Det (B*Bx,, x) < 


Isi,ksan lst, kon 


which, together with (9), yields (7). 
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3. Tneorem 2. Let A, B be two c. c. operators in 9. If td}, ef, 
uit, foi; are the eigenvalues of A*A, B*B,(AB)*(AB), (A + B)*(A + B), 
respectively, then for any two non-negative integers m, n: 


epics SNe tN Kn +] (10) 


/ 
< , 
Vv Im4n+1 > Vv : 1 + V Kn4 1; (11) 


Proof: Let }xi{, }y.{ be two complete orthonormal sets such that 
AA*x, = dx, B*By, = ny, (0 = 1,2,...). Let AB = WI be the polar 
factorization of AB, where // is the non-negative square root of (A B)*(AB) 
and W is partially isometric. For any z of , we have (//z, z)? = (W*ABz, 
z)? < (AA*W2z, Wz)(B*Bz, z). If z is subject to conditions |/z) = 1, 
(z, W*x,) = O(1 <i < m) and (2, vy) = O(1 <j < n), then (AA*W2z, Ws) 
<]]W2i!*\m 4 S Ama (B*Bz, 2) < «x, 4, and therefore (Hz, 2)? < N41. Ket 
Thus we get (10) by applying the minimum-maximum principle to //.. The 
proof of (11) 1s similar. 

The next theorem deals with asymptotic distribution of eigenvalues. 

Tueorem 3. Let A, B be twoc. c. operators in. Let |d{, {mf, fait be 
the eigenvalues of A*A, B*B and (A + B)*(A + B) respectively. If, for 
some exponent r > 0, we have 

lim 2"), = wo < + ~ and lim n'x, = 0, 


n> « n—> © 


then 


lim n’o, = w. 
un x 
Proof: Leth be a fixed positive integer. By (11), we have 
V on4inty S V Nan tj + Vv Kn4 , (" _ Ly (ae te) ao eee) 
which combined with (12), gives 


lim 2’ + 2 "< [(4 + 1) h|’ F.e9'f* (14) 


n—> o 


If we apply (11) to A = (A + B) + (—B), we get inequality V o,, 4; => 
V Nats V ky 41, Which together with (12) yields 


- -/2 , wr /2 ™ 
lim 2”? 9, > [h/(h+1)]'* w (15) 


Then (13) results from (14), (15) by letting h +o. 
4. We turn now toa general inequality for symmetric gauge functions. 
THEOREM 4. Leta; > a. >...>a,2°0,05>h>...>6,>0. Thena 


10 


necessary and sufficient condition that 
Pa, (gy s+ vy ay) < P(h, bo, sey by) (16) 


for all symmetric gauge functions ® of n variables, is that 
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k k 
~as +d (lo k<n). 
s=1 t=] 

The necessity of (17) is immediate, when we consider, for each k, the 
particular symmetric gauge function ® defined by 

P(E, f,..-,&%) = Max  (lé&| + [&| +... +[Eel). 
ii<h< <tk 

For the proof of the sufficiency of (17), we only indicate that it results from 
the following two lemmas. 

Lema 2. Let {,} bea set of n real numbers. LetO < a < land leth, j be 
two distinct fixed indices. I} 


Cn = ab, + (1 — a)d,, c, = (1 — a)h, + aby, c, = b, fort # h, j, 
(1S) 


then ®(¢4, C2, . . «5 Cn) < B(di, be, . . ., bn) for all symmetric gauge functions © of n 
variables. 

Lemma 3. Let {a,}, {b,} be two sets each of n real numbers satisfying (17) 
and a; > a;4,. Then by the successive application of a finite number of substi- 
tutions of form (18), {bi can be reduced toa set {di} such thata,;< dj (1<1i< 
n ) : il 

5. Tueorem 5. Let A, B be two c. c. operators in . If {dj}, {xf and 
{ o,} are the eigenvalues of A*A, B*B and (A + B) * (A + B) respec- 


tively, then 


P( V a1, Wenss V on) < P( Vy V das ’ V vn) + P( Vv Ki, V ke, cee V Ka) 
(19) 


holds for every symmetric gauge function ® of any number n of variables. 
Proof: By Theorem 4, we need only to prove the following special case of 
(19): 


n n n 
T V¥e8 5 Vu tet Va = 1,2,3,...). (20) 
1 t=] ‘ 1 


Consider the polar factorization of A + B = WHI, where // is the non- 
negative square root of (A + B)*(A + B) and W is partially isometric. 
Let 2, 22, ..., 2, be m orthonormal elements such that //z, = V 042; <ts 
n). Wecan find a unitary operator l’ such that (A + B)z, = W*(A 4 
B)z, = Hz, Al<i<sn). Then 


n 


> V4, = dL (U(A + B)z, 2%) = D> (UAx, 2%) + 2 (UBz, 2) 


t= 1 t=] 


and (20) follows from 
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¥ (UAz, 2) 1 < Vy > (UB, %) | < Vk, 
i=] i=] i=l 


t=1 
which are true by (4). 


* This work was supported in part by the Office of Naval Research. 

' Von Neumann, J., “Some Matrix-Inequalities and Metrization of Matric-Space,” 
Tomsk Univ. Rev., 1, 286-300 (1937). 

2 The eigenvalues will always be arranged in order of decreasing absolute value 

* Weyl, H., “Inequalities Between the Two Kinds of Eigenvalues of a Linear Trans- 
formation,” these PROCEEDINGS, 35, 408-411 (1949) 

‘Horn, A., “On the Singular Values of a Product of Completely Continuous Opera 
tors,” [bid., 36, 374-375 (1950). 

® Weyl, H., “Das asymptotische Verteilungsgesetz der Eigenwerte linearer partieller 
Differentialgleichungen,’’ Math. Annalen, 71, 441-479 (1912). 

* Hardy, G. H., Littlewood, J. E., and Polya, G., Inequalities, Cambridge Univer- 
sity Press, Cambridge, 1934, pp. 44-49. 

7 Von Neumann, J., “Uber adjungierte Funktionaloperatoren,” Ann. of Math., 33, 
294-310 (1932) 

* Polya, G., “Remark on Wey!’s Note: Inequalities Between the Two Kinds of 
Eigenvalues of a Linear Transformation,”’ these PROCEEDINGS, 36, 49-51 (1950). 

* Fan, K., “On a Theorem of Weyl concerning Eigenvalues of Linear Transforma- 
tions, II,"’ /hid., 36, 31-35 (1950) 

Following von Neumann, a gauge function & (in the sense of Minkowski) is called 
symmetric, if P(&, &, » én) = Pegs, ert), , €ntj,) for any combination of signs 
e, = +1 and for any permutation (ji, j2, , jn) of (1, 2, ...,m). For general properties 
of symmetric gauge functions, see Schatten, R., A Theory of Cross-Spaces, Princeton 
University Press, Princeton, 1950, pp. 84-92 

'! The proof of Lemma 3 is similar to the corresponding part of the proof of Muir- 
head's theorem (see the reference in footnote 6). Using Lemma 3, one can also prove: 
For two sets fa}, {b;} each of nm non-negative numbers arranged in descending order, 
condition (17) is equivalent to the existence of n? non-negative numbers pi satisfying 


n 
(1) and a = > Pub) Cl <i n). 
j=1 


INVARIANT DIFFERENTIAL FORMS IN SEVERAL GROUP 
VARIABLES AS SOLUTIONS OF PARTIAL DIFFERENTIAL 
EQUATIONS IN FRECHET DIFFERENTIALS* 

By ARISTOTLE D. MICHAL 
CALIFORNIA INSTITUTE OF TECHNOLOGY 
Communicated by G. C. Evans, September 4, 1951 

1. Introduction. There is a particularly interesting class of partial dif- 
ferential equations in Fréchet differentials that characterize abstract in- 
variants and in particular differential forms that are left invariant by the 
parameter groups m and my, of a differentiable abstract group 4. We shall 
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concern ourselves with alternating differential forms in several group 
variables invariant under m, and/or , of an abstract group A imbedded in a 
Banach space B. These groups A are topological groups but they are not 
necessarily compact nor even locally compact. Hence the Haar measure 
theory is not available here in general. Of course if A is assumed to be 
compact or locally compact, then Haar measure theory is applicable and 
leads to useful results in invariant theory. 

The abstract group A will be assumed to be an open set imbedded in a 
Banach space such that the group product a and inverse a’ together with 
the unit «form a group of Class C> (1.e., the group functions a8 and a’ have 
continuous second successive Fréchet differentials throughout A). Other 
assumptions on A will be explicitly stated. A Fréchet differential of a 
function f(a) with arguments and values in Banach spaces will be denoted 
usually by f(a; y) if the increment is y and the differential is ata = a. A 
partial Fréchet differential in the 7th variable of a function f(a, .. ., 
a, ..., a) will be denoted by fa,(ai, ..., ai, ..., Qn; y) if y is the incre- 
ment and the evaluation is at a, = a;. In the special case of a partial Fréchet 
differential in the first variable, we shall often write more simply f(a, (1, ..., 
3,; y) for the Fréchet differential in a at a = a@ with increment y of the 
function f(a, Bi, .. ., Br). 

The parameter groups m and m, are respectively the continuous trans- 
formation groups 


™: a = af (8 is the parameter), 
wT: B = a8 (a is the parameter. 


Let us define two functions 


f(a, a, y) = So’ ay 
def , 


gly, B, o) = yBo 


where y is an arbitrary but fixed element of 4 while a, 8 and o are independ- 
ent variables over A. The two linear differential forms 


(8, 58) = [f(o, B, 7; 5B) leu 
P(a, da) = [g,(y, a, oj ba) |, 24 


(1.1) 


can be shown to be solvable linear functions of their second variables with 
inverses 22> '(8, &) and ®~'(a@, n) respectively. It is to be noticed that the 
notation —1 does not mean group inverse. In fact, the values of the func- 
tions 2~' and #~' are not in general in the group 4. It can be shown that the 
differential forms 2(, 68) and ®(a, da) are invariant under mw. and m re- 
spectively for all 8, a of the group 4. Furthermore, they satisfy the follow- 
ing differential equations in Fréchet differentials' (the generalizations of the 
Maurer-Cartan’ equations of the Lie theory) 
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0( 8, 6:8; 6.8) — 2(B, 62.8; 68) = C[Q(B, 6,8), Q(B, 6.8)] (1.2) 


and a similar equation for (8, 68). The abstract alternating bilinear func- 
tion C(é, 7) is the generalization of the structural constants? of the Lie 
theory. 

Finally, to shorten the statements of our theorems, we need to define one 
more concept and make clear a notation for differential operators. A 
group A will be called d-arcwise connected if A is arewise connected with are 
functions that have continuous derivatives. 

If A(B, &) is a fixed function of Class C; on AB to B, then we define the 
differential expression X;,F as 


X;,F = Fou Bi, Bo, Y, Sy 22, - A( (Bi, &)) + F(~p,, Bo, A( fi, & »y see (1. 3) 


1 
+ F(p,, be, y, A(B., &} 21), 2) + FB, Be, y, 2, A(B2, & 2 


where F( (8), B», y, 21, 22) 1s on A°B* to B,, trilinear in y, 2; and z2, and B, isa 
Banach space, not necessarily distinct from B, while A°B* denotes that 
3), B. range over A and y, 2), 2. over B, 

2. A Typical Theorem on Invariant Two-Point Differential Forms._-The 
following is a representative theorem. We believe that both the results and 
methods of proof are novel even if the Banach spaces B and B, reduce to 
finite dimensional Euclidean number spaces (the same can be said for 
Theorem I] of §3). Only a few remarks on the method of proof can be given 
here as the complete proof is long. The results are global, i.e., in the large. 

THeoreM I. Let A ¢ B,a Banach space, be an open d-arcwise connected 
group of class Cy. A necessary and sufficient condition that a two-point dif- 
ferential form W(B,, Be, 681, 68», 628») of class C, (for By, Bo in A, and 6B, 6; B» 
in B) with values in a Banach space B,, not necessarily distinct from B, linear 
in 63, and alternating bilinear in 6,82, 6282, be invariant under the second 
parameter group m. of A is that W satisfy the following partial differential 
equation in Fréchet differentials 

X,W = 0 (2.1) 
with the understanding that the &~'( 8, £) defined in connection with (1.1) 1s 
the fixed function A(B, £) in XW. 

For the sufficiency part of the proof, it is not difficult to establish first 
that there exists a function A(B), 8», £1, m, m2) on A®B* to By, linear in &, and 
alternating bilinear in y and 7» such that 


W (Bi, Be, 681, 618, 6282) = ACP, Be, (Bi, 681), LC Be, 182), Q(B», 6282)). 


(2.2) 
The theory of the abstract Pfaffian differential equation’ 


P(6, &; 68) = C(P(8, &), 2(B, 6B)) 
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satisfied by P(g, £) = 2(8, ®-'( 6, &)) for all 8 « A and £, 638 « B, and of the 
abstract® equation (1.2) are used to show that A( 6), 6», £m, m2) must be in- 
variant under m as a function of the group variables 6, and 6, whenever 
(2.2) satisfies the equation (2.1). Finally it is shown that 


W = X( Bi’ Bo, {2( B, 62B)), q2( Bo, 6; Bs), (2( Bo, 6B») ) (2. 1) 


for some function A(B, & m, m2) of only four variables. From this the suf- 
ficiency part of the condition follows without much difficulty. The proof 
for the necessity part of the condition uses a method that essentially re 
verses the steps for the sufficiency proof but does not use the d-arewise con- 
nectedness of the group A. 

There is a corresponding theorem for invariance under m with the 
corresponding partial differential equation 


(2.5) 


in which A(, &) is now replaced by 2-'({, &). 

3. A Typical Theorem on Invariance Under Both Parameter Groups. 
Necessary and sufficient conditions that a differential form HW’ be invariant 
under both m, and m2 can clearly be given in terms of HW’ satisfying the pair of 
partial differential equations (2.1) and (2.5). However, there is another 
type of theorem proved by the writer in which the necessary and sufficient 
condition for invariance is not in terms of solutions of partial differential 


equations. For example, a necessary and suthicient condition (under the 


, 


restrictions of Theorem 1) that W be invariant under both rm, and mo ts that VW 
be of form (2.4), that (6, & m, m2) be a class function of the group variable B 
satisfy the functional equation 


d(B, C(é1, £), m, m2) + ACB, £1, Cm, &), m2) + ACB, &, m, Cm, &)) = O 
(3.1) 


forall B «e Sand &, &, m, ne € B. 
There is still another type of theorem that combines some of the features 
of the last theorem together with those of Theorem I. We shall state it as 
THEOREM II. Under the restrictions of Theorem I, a necessary and suf- 
ficient condition that the two-point diffierential form W be invariant under both 
parameter groups mw, and m2 of the group A is that W satisfy the abstract partial 
differential equation (2.1) and that the function defined by 


MB, 1, m, m) = Wie, B, Qe, &), 2-8, m), 2-(B, m)] 


(e, the unit of A) be a class function of B and satisfy the functional equation 
(3.1). 

Now, in the theory of abstract one-point alternating differential forms 
invariant under both m, and wz (in the interim of publication’), the abstract' 
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Poincaré-Stokes differential of the form is the zero element of the value 
Banach space B,. It can be shown by counter examples that such is nol 
in general the case for mu/tiple-point forms invariant under both x and 7, 
i.e., the abstract Poincaré-Stokes differential may not be zero with respect 
to any of the group points. 

4. Concluding Remarks.--We mention that the theory of analytical 
functions of abstract®*® variables in complex Banach spaces (as well as real 
Banach spaces) has immediate application in case we take the spaces B and 
B, to the complex Banach spaces. One can also consider invariant mul- 
tiple-point Hermitean differential forms. 

The author's differential calculus’ for functions of topological abelian 
group variables can be used in generalizing the subject of abstract differenti- 
able groups and their invariant differental forms whenever general topologi- 
cal abelian groups are taken in the place of the Banach spaces. Some theo- 
rems that involve d-arewise connectedness would need modification. If 
the group A is imbedded in a locally convex linear topological space, then 
the theorems of this note continue to hold with the understanding that the 
notion of a function of Class C, is made a little more restrictive. 

In conclusion, we wish to mention that these studies on invariants have 
connections with harmonic forms and their generalizations, classical as well 
as modern differential geometry, and the author's general differential 


geometries with abstract coordinates.* 


* This note represents a summary of a chapter of an extensive monograph in prepara - 
tion: Michal, A. D., Differential Equations in Abstract Spaces with Applications. We 
shall refer to this by the abbreviation DEASA. This note can also be considered as a 
continuation of the author's early studies on non-additive integral invariants: Michal, 
A. D., Functionals of R-Dimensional Manifolds Admitting Continuous Groups of Point 
Transformations, Trans. Am. Math. Soc., 29, 612-646 (1927). 

' Michal, A. D., and Elconin, V., Differential Properties of Abstract Transformation 
Groups with Abstract Parameters, Am. J. Math., 59, 129-143 (1937); Michal, A. D., 
Differentiable Infinite Continuous Groups in Abstract Spaces, Revista De Cienc: 
(Lima, Peru), 50, 131-140 (1948); Michal, A. D., DEASA, Chapter 1X 

* Cartan, Elie, La Théorie Des Groupes Finis Et Continus Et La Géométrie Differentielle, 
Gauthier-Villars, Paris, 1937; Eisenhart, L. P., Continuous Groups of Transformations, 
Princeton Univ. Press, 1933 

§ Michal, A. D., DEASA, Chapter IX. 

* Michal, A. D., General Tensor Analysis, Bull. Am. Math. Soc., 43, 8394-401 (1937). 

® Michal, A. D., and Martin, R. S., Some Expansions in Vector Space, J. Math. Pures 
Appl., 13, 69-91 (1934) 

® Michal, A. D., DEASA, especially Chapter I] and Chapter III (some of the un- 
published 1931-1932 work of Michal-Martin is included here); Michal, A. D., On a 
Non-Linear Total Differential Equation in Normed Linear Spaces, Acta Mathematica, 
80, 1-21 (1948) 

’ Michal, A. D., Differentials of Functions with Arguments and Values in Topological 
Abelian Groups, these PROCEEDINGS, 26, 356-359 (1940); Michal, A. D., Functional 
Analysis in Topological Group Spaces, Math. Magazine, 21, 80-90 (1947) 
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* Michal, A. D., General Differential Geometries and Related Topics, Bull. Am. Math. 
Soc., 45, 529-563 (1939); Recent General Trends in Mathematics, Science, 92, 563-566 
(1940); Infinite Dimensional Differential Metrics with Constant Curvature, these 
PROCEEDINGS, 34, 17-21 (1948). 


REPRESENTATION OF NUCLEON-NUCLEON SINGLET S 
SCATTERING* 
By G. Breit AND M. C. Yovits 
YALE Universiryt 
Communicated September 27, 1951 


In early work! on proton-proton scattering it has been found convenient 
to deal with the function 
f = {2e/[—1 + exp (2xn)|} cot Ay + 4y — 2 + 2R.P.[T’(in) /T(in)} 
—2Inn 


where 


= Euler-Mascheroni constant, 


Ky is the phase shift, Z,Z’ are atomic numbers of the colliding particles and 
I(x) is the gamma function of x. The special reason for introducing the 
function f was the fact that the logarithmic derivative of rR = F, where R 
= radial function, at the outer boundary of a short range potential well is 
expressible in terms of f and quantities which are power series in the energy 
I. For potentials of very short range the function / is' a constant and the 
slope of f plotted against / is proportional to the range of force as is clear' 
from the formulas connecting / with the logarithmic derivative and other 
quantities. While the function / and its approximate linearity have been 
used as intermediate steps in the earlier calculations a strong impetus to its 
adoption in recent years has been given by Schwinger? who found it con 
venient to deal with this quantity by a process of conversion of the radial 
differential equation to an integral equation and representing f as a varia- 
tional integral. Following Schwinger, Blatt? and Jackson and Blatt*® have 
employed / and its limiting form for ZZ’e? = Oas a convenient means of ana- 
lyzing experimental material while Breit and Bouricius* have made similar 
studies without the variational principle. Other treatments of the simpler 
properties of / have appeared* since; in these the variational method is 
avoided and the first two terms of / are dealt with directly by means of the 
differential equation. The impression has nevertheless been created that 
the variational method puts f in its true (‘natural’) light and that its 


power series representation can be studied by means of the variational 
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method more simply than otherwise. It is apparent, on the other hand, 
that in a problem which amounts to a study of an ordinary differential 
equation of second order all methods must have rather immediate connec- 
tions. It is the modest purpose of the present note to describe an approach 
in which the interrelations of the variational and direct treatments can be 
seen in a simple manner and to point out a combined use of the two ways if 
one wishes to obtain a simple expression for the general term in the ex- 
pansion of f in powers of /: or else in powers of a linear parameter such 
as a change in depth of a potential energy well. 

It follows from the stationary property of the variational integral that an 
error in the wave function of relative order ¢ produces an error in f of rela- 
tive order e or higher. This fact is readily seen to indicate the existence of 
identities between integrals containing coefficients of // in the expansion of 
F in powers of £. The argument from here on is essentially one regarding 
the form of the answer rather than the degree of approximation. It is con- 
venient to start out with the slight generalization of Schwinger’s form for 
0//O0F which has been used by Breit and Hatcher,‘ viz. 

of — f’)/ox =a f-" (u’? — u?) dr (1) 


where / is the radius of the potential well, /’ is f for a comparison potential, 
x is k® where k is the usual 2x divided by wave-length, a = h?/ye?ZZ’ is the 
Bohr length with the reduced mass denoted by yw and u, u’ are suitably re- 
normalized functions F for the actual and the comparison potentials. The 
normalization of u, uw’ is such that if the function is continued as the solution 
of the Coulomb wave equation to r = 0 the value of the function is unity. 
The form of equation (1) is identical with that of Schwinger but it is used 
here for a variable 4 0. One expands 


u = >> (x — K)"U, (1.1) 
0 


and obtains from (1) 
ad (x = Ko)" do SK? (uu, — ua; dr = 


cd 
rs n(x — Ko)®—"[{ f™ 


” ee | n 
n 


f ™ es (x = Ko) "f{™, 


and a similar formula for /’ define the expansion coefficients [(", f". It 


follows from equation (1.2) that 


nf f~™ —f'™) =ad> SL” (u'u', — um, dr. 
! J a i | =F 


i+j=un 1 
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On the right side of this formula there enter all the expansion coefficients 
Uy, Uy, . . «» Uy. The variational principle shows, however, that if n = 
2v + 1 is odd then the knowledge of um, . . ., u, suffices for the determination 
of f,..., ft) = f™ so that the entrance of u,,),... 4,1 in equation (2) 
is unnecessary. Similarly if m = 2viseven then the knowledge of uw, ... 4, 
should suffice and the entrance of u,,;, .. . M,—1 iS unnecessary. There must 
exist, therefore, identities between the integrals in equation (2). From here 
on one does not need the variational principle because one finds from the 
differential equation that 


(s— 1,#+ 1) = (s, 2), (2.1) 
where 
(st) = SO" (u’, wu’, — u,u,) dr. 
It is thus seen that 
fm — f° = a(s, s), (2 = 2: 
fm —f'™ = a(s — 1, s), (n =: (2.3) 


which is a general answer to the problem of obtaining the coefficients in the 
expansion of f in powers of /:. For a linear parameter \ entering the dif- 
ferential equation as in 

d?u/dr* + («x + -.../% = (23) 


one obtains 


of/OX = —a a gu? dr. 


The introduction of the symbol 
Is, th = A” tattag dr 
yields, by an argument very similar to that for « — ko, 


f = fy + fin + A°foy + 


-alfs — 1, s], (wn = 2s). (4.1) 


The above discussion applies to other problems than the nucleon-nucleon 
singlet S scattering although the more important applications of the f func 
tion have been so far principally in this field. It is seen that the technique 
of calculation is simpler by the direct approach but that the stationary 
property of the variational integral, while not indispensable, is nevertheless 
a valuable guide. 

* This paper has been read at the Washington Meeting of the American Physical 


Society, Spring, 1951. 
t Assisted by the joint program of the ONR and AEC 
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